Bosonic effective action for interacting fermions 
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We compare different versions of a bosonic description for systems of interacting fermions, with 
particular emphasis on the free energy functional. The bosonic effective action makes the issue 
of symmetries particularly transparent and we present for the Hubbard model an exact mapping 
between repulsive and attractive interactions. A systematic expansion for the bosonic effective 
action starts with a solution to the lowest order Schwinger-Dyson or gap equation. We propose 
a two particle irreducible formulation of an exact functional renormalization group equation for 
computations beyond leading order. On this basis we suggest a renormalized gap equation. This 
approach is compared with functional renormalization in a partially bosonized setting. 



PACS numbers: ll.10.-z, 11. 10. Hi, 11. 10. St 
I. INTRODUCTION 

Systems of strongly correlated electrons often show the 
phenomenon of spontaneous symmetry breaking or a be- 
havior close to it. Examples are antiferromagnetism or 
superconductivity in the Hubbard model [1]. Sponta- 
neous symmetry breaking (SSB) is also common to many 
other fermionic systems, like effective quark-models for 
the strong interactions [2] or ultracold fermionic atoms 
[3] . The spontaneous breaking of symmetry is introduced 
by a bosonic order parameter or vacuum expectation 
value (vev), typically corresponding to the condensation 
of fermion-fermion or fermion-antifermion (electron- hole) 
pairs. Examples are Cooper pairs ~ ipip for supercon- 
ductors or the spin vector rh ~ iprip for ferromagnets. A 
reliable quantitative description of SSB - for example the 
computation of the order parameter and the mass gap 
in the low temperature phase - is not easy in a purely 
fermionic language. The simplest possible effective po- 
tential for the magnetization in a ferromagnet involves 
already four powers of m, e.g. 



U = -ym + - [m ) 



(1) 



The term ~ m corresponds to an eight-fermion- 
interaction which is hard to compute in a purely 
fermionic picture. Moreover, the effective potential (or 
free energy) often does not take a simple polynomial form 
as in the example (1). 

Standard approaches to this problem are the mean field 
theory (MFT) or the solution of a gap equation based on 
the Schwinger-Dyson (SD) equation [4]. In both meth- 
ods one treats effectively the fermionic fluctuations in a 
bosonic "background field" or "mean field" . This is man- 
ifest in MFT where the fermionic fluctuations are trun- 
cated in quadratic order. The remaining fermionic func- 
tional integral is Gaussian and the free energy for a given 
mean field can be computed in this approximation. Mini- 
mization of the free energy yields a self-consistency equa- 
tion for the mean field. The shortcoming is the complete 
omission of the effective bosonic fluctuations or, equiv- 
alently, the fermionic fluctuations beyond the quadratic 



approximation. For strongly interacting systems this ap- 
proximation is not expected to be quantitatively accu- 
rate. 

This issue is most easily addressed by partial bosoniza- 
tion via a Hubbard-Stratonovich- transformation [5]. In 
the partially bosonized language the partition function 
is written as a functional integral over fermion and 
(composite) boson fields and MFT precisely corresponds 
to the neglection of the bosonic fluctuations. Partial 
bosonization is, however, not unique and the results on 
SSB depend strongly on the choice of the mean field. 
This is a well known problem for the Hubbard model 
where certain implementations of "partial bosonization" 
are consistent with the spin rotation symmetry but fail 
to reproduce the Hartree-Fock result (i.e. lowest order 
Schwinger-Dyson equation) whereas others are consis- 
tent with Hartree-Fock but do not exhibit the correct 
symmetries . The problem can be traced back to the 
possibility of Fierz-reordering of a local interaction. For 
the Hubbard model a systematic discussion of the strong 
quantitative importance of the "Fierz ambiguity" can be 
found in [7], In a more general context the origin and 
the cure of this "Fierz ambiguity" are discussed in detail 
in [8]. There it is shown that even a rough inclusion of 
the bosonic fluctuations greatly reduces the ambiguity. 
The remaining ambiguity can be used for an error esti- 
mate of a given approximation - physical results must, of 
course, be independent of the choice of a specific partial 
bosonization. All this clearly demonstrates the impor- 
tance of the bosonic fluctuations neglected in MFT. 

The SD-approach is formally a fermionic formulation 
and the lowest order gap equation is one loop exact. 2 
However, an apparent problem is the reconstruction of 
the free energy from the solution of the gap equation. 
This becomes crucial when the gap equation admits so- 
lutions with different order parameters and the free en- 
ergy corresponding to the different solutions has to be 



1 This situation was summarized by 
based on nonlinear fields is proposed. 
2 This contrasts with MFT [8]. 



where also a cure 
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compared. 3 Furthermore, the lowest order SD-equation 
for a given order parameter gcnerically takes the form 
of a MFT computation for a particular choice of partial 
bosonization. This shows that the gap equation leaves 
out the effective bosonic fluctuations 4 , just as MFT. This 
shortcoming becomes particularly important when the 
composite bosons are correlated on large length scales 
(or have a small renormalized mass) , as characteristic for 
critical behavior in the vicinity of a second order phase 
transition. 5 

Going beyond the MFT and SD approaches, which are 
established since many years, requires a method which 
should (a) allow for a computation of the free energy 
(b) be one loop exact and (c) include the effects of the 
bosonic fluctuations. One possible approach relies on the 
exact renormalization group equation for the effective av- 
erage action [10]. It is based on partial bosonization, but 
the Fierz ambiguity of MFT can now be cured [8] by 
the inclusion of the bosonic fluctuations and an appro- 
priate "rebosonization" of multi-fermion interactions at 
every scale [11]. Partial bosonization brings a redun- 
dancy of the description: a four fcrmion interaction am- 
plitude can be composed from a direct fermionic vertex 
and a piece from the exchange of bosons. Rebosoniza- 
tion avoids this redundancy by eliminating (parts of) the 
direct vertex in terms of interactions involving bosons. 
Partially bosonized functional renormalization has been 
employed successfully [12] to the issue of chiral symme- 
try breaking within effective interacting quark models for 
the strong interactions. Recently, it has also led to a de- 
scription of the low temperature antiferromagnetic phase 
for the two dimensional Hubbard model for small doping 
[13]. 

In this paper we compare partial bosonization with a 
method based on the "bosonic effective action" (BEA). 
The latter approach was pioneered in solid state physics 
by Luttinger and Ward [14]. (See [15] for a review.) As an 
alternative to partial bosonization one introduces bosonic 
sources for fermion bilinears, but no explicit bosonic fluc- 
tuating fields. As usual, the effective action obtains from 
an appropriate Legendre transform and is a functional 
of bosonic fields. We argue that there is no need to in- 
clude fermionic sources - as a result the BEA depends 
only on bosonic variables. The original fermionic prob- 
lem is then completely translated into a purely bosonic 
language. In particular, this avoids the complications of 
redundancy by construction. Indeed, all operators with 
possible nonzero expectation values are of bosonic nature 



3 See [9] for a recent discussion of this issue in color super- 
conductivity. Their attempt to reconstruct the free energy 
may not always work. 

4 In contrast to MFT the bosonic fluctuations do not con- 
tribute in lowest order in the coupling constant. 

5 The gap equation only gives the MFT values for critical 
exponents. 



- as, for example, the fermion propagator or any order 
parameter. The general possibility of a purely bosonic 
description should therefore be of no surprise. It is also 
clear that a purely bosonic description will, in general, 
not only involve local fields. The bosonic fields corre- 
spond to fermion bilinears V«V'/3 where ip a and ipp may 
be fermion fields at different locations. 

The BEA is based on the two particle irreducible (2PI) 
effective action [14,16,17]. Indeed, the loop expansion 
of the BEA involves 2PI graphs in the fermionic lan- 
guage. However, we choose here a complete translation 
to a bosonic picture, such that the BEA depends only 
on composite bosonic fields. Some of the formulae in the 
first sections of this paper are well known from earlier 
work on the 2P7-effective action. The present derivation 
highlights the close link between the S'D-equation and 
the bosonic effective action and makes several aspects 
particularly simple and transparent. 

The new aspects of this work concern several issues. 
We put emphasis on the free energy functional which is 
needed for a comparison of different solutions of the gap 
equation. In particular, we present an explicit local for- 
mula in the case of a local gap in sect. VII. This can 
be generalized to multi-fcrmion interactions (beyond four 
fermion interactions) in a straightforward way [18]. (Pre- 
vious discussions of the effective potential or free energy 
can be found in [19]). We also construct a systematic 
mapping between the fermionic effective action (IP/) 
and the bosonic effective action (2PI). This helps for 
an understanding of the issue of redundancy. In partic- 
ular, an exact S-D-identity for the BEA permits simple 
systematic expansions. 

Furthermore, we provide for a general discussion of the 
symmetries of the BEA. They go far beyond the sym- 
metries of the fermionic effective action since the BEA 
respects all symmetries of the interaction term. In par- 
ticular, two models with the same interaction but dif- 
ferent fermion kinetic or mass terms lead to the same 
BEA. Their difference only results in a different source 
term in the field equations. We show how the symme- 
try transformations can be used to map models with dif- 
ferent fermionic quadratic terms into each other. The 
large degree of symmetry is particularly apparent in a for- 
mulation where particles and antiparticles (electrons and 
holes) are treated as different components of a fermionic 
Grassmann field. This formulation leads also to simple 
mappings between models with different interactions. As 
an example, we discuss an exact mapping between the 
Hubbard models with repulsive and attractive interac- 
tion. Results from numerical simulations in the attrac- 
tive case can therefore be used as an information about 
appropriate quantities in the model with a repulsive in- 
teraction. 

One of the main tasks of this paper is a detailed com- 
parison between the BEA and partial bosonization based 
on the Hubbard-Stratonovich-transformation. We con- 
centrate on the correlation functions that are computed 
in both methods. This highlights the advantages and 
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problematic points of the two approaches. 

For a first time, we derive exact functional renormal- 
ization group equations within the 2PI formalism. We 
restrict our aim here to a formal establishment of these 
equations and a first qualitative discussion, while post- 
poning practical applications to future work. In partic- 
ular, we propose a simple renormalized gap equation. It 
has the same structure as the lowest order gap equa- 
tion, but with "microscopic" couplings and propagators 
now replaced by renormalized ones. Even for local mi- 
croscopic interactions the renormalized coupling will not 
remain local, in general. This can lead to new structures 
in the space of possible solutions. For example, one may 
find a gap equation for a nonlocal superconducting order 
parameter in the Hubbard model which is not present 
in lowest order. We also propose an improvement of the 
"lowest order" flow equation for the four fermion ver- 
tex by including the effects of a gap. This can avoid 
the previously observed divergence of the renormalized 
coupling in the Hubbard model [20-22]. Furthermore, 
we develop "bosonic renormalization group equations" 
which permit to establish a close contact with the renor- 
malization group for bosonic systems. This allows for a 
description of critical exponents at second order phase 
transitions beyond the lowest order "mean field values" . 
We compare the different versions of the functional renor- 
malization group for the BEA with the one for partial 
bosonization. 

This paper is organized as follows: After defining the 
BEA in sect. II we translate in sect. Ill the fermionic 
SD-equation into a bosonic formulation. The resulting 
identity for the BEA can be used for the construction of 
a systematic loop expansion for the BEA and the free 
energy. The final results of these sections are well-known 
[15]. It is obvious that the BEA is a rather simple object 
which can explicitely be written down in rather high or- 
ders of the loop expansion. The price to pay is the com- 
paratively high complexity of the field equation which 
needs to be solved in order to find the minima of the free 
energy. The BEA also has a very high degree of sym- 
metry (sect. IV). It is invariant under all anomaly free 
transformations leaving the "classical" fermionic inter- 
action invariant, irrespective of the form of the classical 
fermion propagator. 

In sects. V-VIII we turn to models with local four- 
fcrmion- interactions. Based on the well-known local gap 
equation (sect. VI) and a simple local formula for the 
free energy (sect. VII) we discuss antiferromagnetism in 
the two dimensional Hubbard model in sect. VIII. Al- 
ready the lowest order computation of the free energy 
realizes all symmetries in a simple linear way (includ- 
ing spin rotations) and is nevertheless consistent with 
Hartree-Fock. We compute the free energy for the anti- 
ferromagnetic state and sketch how it can be compared 
with other phases. 

In sects. IX, X we turn to the symmetries of the Hub- 
bard model and describe the exact map between the mod- 
els with repulsive and attractive interaction. 



Sect. XI is devoted to partial bosonization. We first 
describe a mixed effective action which supplements the 
BEA by a redundant but essentially local piece for the 
propagation and interaction of fermions. Subsequently, 
we turn to a "local fermion boson model" as defined by 
a Hubbard-Stratonovich transformation. The identities 
relating it to the mixed effective action or to the BEA 
are exact. However, the different versions of bosonization 
will typically use different approximation schemes. 

The powerful tool of the exact renormalization group 
equation for the effective average action [10] can be di- 
rectly implemented for the BEA (sect. XII). There is 
no need for rcbosonization [11] in this formulation. In 
particular, we advocate the use of a "renormalized gap 
equation" in sect. XIII and propose a "gap-improved" 
flow equation for the evolution of an appropriate running 
"four fermion coupling" . We derive the bosonic func- 
tional renormalization group equation in sect. XV. After 
comparison with the functional renormalization group for 
partial bosonization in XVI our conclusions are drawn in 
sect. XVII. 



II. BOSONIC EFFECTIVE ACTION 



We start by collecting all fermionic degrees of freedom 
in a set of Grassmann variables ip a , where the collective 
index a labels momenta (or locations) as well as possible 
internal degrees of freedom including spin. It also dif- 
ferentiates between what is usually called tp and ip, e.g. 
electrons and holes or particles and antiparticles. For 
a model with quartic fermionic interaction the partition 
function reads 



'hj} = J 



Dip exp(rj a ip a + ^j a p4> a 4>p 



"24^ a/J7 /37 



(2) 



Due to the anticommutation property of the Grassmann 
variables A a( g 7< 5 is totally antisymmetric in all indices. We 
treat here the quadratic terms in the fermionic action as 
a bosonic source term j a p = —jpa and consider first Z as 
a functional of arbitrary fermionic and bosonic sources rj 
and j. With W[r/, j] = In Z[r], j] and using 



{ipa) =ipa = 



dW 



(3) 



we can write the fermion propagator as 



,~ ~ . dW d 2 W dWdW 

Thus W obeys a nonlinear functional differential equa- 
tion relating the dependence on j to the one on 77. Simi- 
larly, we can write 



3 



d 2 W 



dja/3 dj 7 s 



+ ■ 



d 2 W d 2 W 



with 



(2) 



drj a drjp d^li dm drjp drj 7 drj a drjs 
drj a dr] 7 dr)pdi] 5 



d 3 W dW d 3 W dW 



(5) 



drjadrjfjdrij drjs drj a dr]pdris dr) 7 

d 3 W dW d 3 W dW 
+ 



drjpdri-fdris dr\ a dr} a dr] 7 dris drjp 
d 2 W dWdW d 2 W dW dW 



+ 



drjadr]^ drjp drjs drj a dr]s dr\p dr] 7 
d 2 W dWdW d 2 W dW dW 
drjpdrjj dr\ a dr/s dr/pdr]s dr\ a dr/^ 



(6) 



vanishing for rj = 0. 

The fermionic effective action in presence of bosonic 
sources j obtains by a Legendre transform with respect 
to the fermionic variables 



r*#,j] = -W[ri,j] + Va*Pa 
It obeys (at fixed j) 
dT F / d 2 T F 



(7) 



with 



\ /at 



(2) 



dip a dipp 



d 2 W 



(wpy 1 (8) 



a p dr\ a dr}p ' 



(9) 



We may express IV as a functional integral in a fermionic 
background field ip, 



T F [iP,j] = - In / Ehpexpi-Sjty + V] + r)a*Pa) (10) 



where we have defined 

Sj[ip] = ~j a pi> a il>/3 + -^Kp-ySipaippip'yips (11) 



The bosonic effective action for the composite field G a p 
is constructed similarly. For the purely bosonic formula- 
tion motivated in the introduction we restrict the formu- 
lation to vanishing fermionic sources i] a = 0. The BE A 
is then defined as 

r B [G] = -W[O,j] + ^j a pG a0 (12) 



and obeys 



dG 



ap 



(if) 

\ /a 



d 2 r, 



/3,7<5 dGafjdGyS 



a/3,7<5 



where 



= d 2 w[o,j] 

ap tl 8 OjapdjjS 



(13) 



(14) 



We note that I> is an even functional of ip and contains 
precisely the same information as Tb- The bosonic and 
fermionic effective action are connected by relations of 
the type (4), (5). Evaluating them at ip = relates G a p 
to the fermion propagator for vanishing ip, 



G a p = (rf^- ] 



F >ap 



(15) 



and connects the full bosonic propagator to appropriate 
(functional) derivatives of r> 



(2) 



a/3, 7(5 



dlp a > dip p' dipy dtps* 



/ r (2)\-i (r^h- 1 (r^h- 1 (r^h- 1 

\ F ) aa' V F )pP'\ L F J 77 'l i F ) 55' 



i_ I T -1 (2) \ 1 / T"i (2) \ 1 rr (2),_i rr (2),_l 
~r\ L F )p~/\ 1 -F )a8 V 1 F )aj\ L F ) PS 



(16) 



The translation rules (15) (16) - and similar rules for 
higher vertices - permit a mapping between the fermionic 
effective action I> and the bosonic effective action Tb- 
We emphasize that G a p carries two fermionic indices and 
is therefore, in general, a bilocal bosonic object. 

The BEA is simply related to the Helmholtz free en- 
ergy F 



T B = - 



(17) 



The grand canonical partition function Z(T, fi,V) de- 
pends on temperature T, volume V and the chemical 
potential ri which is a part of the source j a p. Similarly, 
the particle number density is a particular combination 
of G a p. It is sometimes convenient to introduce an aux- 
iliary thermodynamic potential 



1 F 

T A G 'i] = T B - ^j a pG a p = — 



Ui 



T7 T J ( 18 ) 



6 Note that due to the antisymmetric index notation 
the bosonic matrix multiplication reads (F ■ G) a p :a ipi = 
^F al 3 n sG jS , a 'i3' with unit operator I a p ta ipr = <W<W - 



2 ' 

5 a0 i5 0a i. This defines {Wg')' 1 . For clarity we often de- 
note contraction of bosonic indices with a dot in order to 
distinguish it from fermionic index contraction, i.e. j ■ G = 
\upG a ? = -iTr(jG). 
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such that the equilibrium states correspond to the ex- 
trema of Tj. We will call Tj a "free energy" as well, 
Tj = F/T. The pressure p, the energy, entropy and 
particle number densities e, s and n can directly be 
gained from the minimum of Tj, i.e. Tj min — —W = 



{F-»N)/T,U. 



U as 



s = 



-U , e = U 
dU 



au du 

dU 



(19) 



Of course, Tb and Tj do not only depend on particle 
number. They are functionals of all possible fermion 
propagators G a p, whereby the physical propagator is 
given by the minimum of Tj . 



III. SCHWINGER-DYSON EQUATION AND 
LOOP EXPANSION 

In this section we translate the Schwinger-Dyson equa- 
tion from the fcrmionic formulation into a corresponding 
exact identity for Tb- The lowest order approximation to 
the Schwinger-Dyson equation results in a closed equa- 
tion relating G to j. This allows us to construct ex- 
plicitely Tb and the free energy. An iterative solution of 
the SD-equation corresponds to a loop expansion for T b ■ 

Differentiating eq. (10) with respect to ip yields the ex- 
act Schwinger-Dyson equation for the fermionic effective 
action 

QT p 1 

= jaf31pa ~ g Kf3~/5(i> a 1p-y1p5) 



-(r^)- 4 V 7 + (r^)^Va 
-ir^V 1 r r (2)\-i (r (2) r i 

\ F ) act 1 \ F )j-y'\ L F >65> > 



d 3 Tj 



dipa'dipy'dtps 



-}. (20) 



The fcrmionic Schwinger-Dyson equation relevant for our 
purpose obtains by evaluating a further differentiation 
with respect to ip at V = 



(lf>) 



1 



a/3 



Z \ / jS 



d 4 T F 



dipadipedtpydips 



S'S 

(21) 



This identity can directly be translated into the 
bosonic formulation by use of the identities (15) and (16): 



Y, 



a/3 



z —ja/3 + -^^afiisG^g + Y a p 

-\ v PjsG a } (Tb^c^jS + G^G v s - G e sG vl 



(22) 



The exact identity (22) is a central equation for this pa- 
per. Its iterative solution will lead to a loop expansion 
for the BEA. All higher order loop effects are contained 
in the term Y a p. An estimate of its size can therefore be 
used as a error estimate for the lowest order approxima- 
tion. 

A first nontrivial approximation neglects the term Y a p 
which corresponds to the last term in the fermionic 
Schwinger-Dyson equation and is of the order A 2 . For 
given j and A the lowest order Schwinger-Dyson equa- 
tion is a closed equation for G 



G a p — —ja/3 + -^^afljsGjS 



(23) 



The solutions of eq. (23) correspond to extrema of the 
free energy F in the lowest order approximation. With 



Ja/3 



dG 



a0 



we can explicitcly compute in the lowest 



order approximation by integrating the differential equa- 
tion 



or, 

8G a R 



^a/3-fSGjS 



(24) 



One obtains (note dG a p/ 'dG^s — Sa-ySps — daS<ipy) the 
two loop effective action 



ri = -Tr In G + -\ aPl sG a pG lS + c (25) 

For qualitative estimates of the properties of a 
fcrmionic systems the lowest order BEA (25) will often be 
sufficient. It is instructive to recall the respective role of 
the two terms on the r.h.s. of eq. (25). The second term 
corresponds to the classical action (without the source 
term ~ j), whereas the first one reflects the fermionic 
fluctuation determinant (^Tr In G = — ^ In dct G _1 ). 
In a graphical representation the first term corresponds 
to a closed fermion loop, whereas the second term in- 
volves two loops connected by the vertex A. Here every 
factor G a p corresponds to a fermion line. A differenti- 
ation of with respect to G corresponds to cutting a 
fermion line and lowers the loop order. The two loop 
BEA is therefore equivalent to the one loop SD-equation 
for the fermion propagator. It may be surprising that 
the classical contribution appears formally as a two loop 
expression. We will see in sect. VII that actually no 
explicit loop calculation is needed for its evaluation. 

It is straightforward to compute the value of the ther- 
modynamic potential F — /iT or U for given values of the 
source j. Inserting the solution (23) the effective action 
can be expressed as 

Ti = ^Tr In G + \j a0 G a p + c! (26) 

where G should be interpreted as a functional of j. This 
yields the lowest order formula for the free energy 
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F-HN 1 

j, - 1 1 - -^a(i3a(i 



-Tr In G - -G a pj af} + c' 



-Tr In G - -\ a p 1 sG a pG 1 s + c 



(27) 



If the field equation (23) admits more than one solution 
the expression (27) can be used in order to determine the 
one which corresponds to the lowest value of the thermo- 
dynamic potential F — /iT. Our formalism therefore is 
well suited to describe competing order parameters and 
first order phase transitions. 

A systematic loop expansion for the BEA can be con- 
structed by an iterative solution for Y a p . In the lowest 

order approximation one finds from differentiation 
of cq. (25). This yields 



(if)) 



a/3, 7(5 



S~1 — 1 1 

~ <J c*7 <J >5 



G a5 G p 1 + A "/37<5 ( 28 ) 



and therefore 



,(2) 



(29) 



Inserting eq. (29) into the formula for Y a p in eq. (22) 
one obtains the two loop expression which corresponds 
to the "setting" sun diagram 

Y a p = —^arpa^flrjjsGrrjGp^GcrS- (30) 

The three loop effective action therefore receives addi- 
tional quartic bosonic interactions 

1^2 = Tl — -^^arpcr^flrijsGafjGrffGpjGcrS (31) 

The approximation (31) is valid as long as |AG 2 <C 1 (in 
a parametric sense). In the opposite extreme |AG 2 | » 1 
we may use a "strong coupling expansion" which will be 
discussed elsewhere. 

Let us finally define the two particle irreducible (2PI) 
part 



f = T B - -Tr In G 
2 



(32) 



It consists precisely of the sum of 2PI diagrams on the 
fermionic level. This is linked in a very direct way to 
the observation that the self energy in the SD equation 
(21) for fermions is one particle irreducible. Differentia- 
tion with respect to G cuts a line and therefore reduces 
the degree of reducibility precisely by one. The one par- 
ticle irreducibility of the self energy dT/dG = T^) +j 
is well known. This would be a contradiction unless T is 
2PI. The proof of two particle irreducibility of V becomes 
therefore extremely simple in our formulation. 



IV. SYMMETRIES 

Many properties of T[G; A] can be understood in terms 
of symmetry transformations. Let us consider a general 
linear transformation (with complex t a p) 



(33) 



We will assume that t a p is regular such that the inverse 
t~p exists. If we also transform the sources j and the 
quartic coupling A according to 

jaP = (t~ 1 )a> a (t~ 1 )p>pj' alf 3 l 
Kp~/S = (t~ 1 )a>a(t~ 1 )p>p(t~ 1 )y~ / (t- 1 )s>sX a >p>ys> (34) 



the action Sj ((eq.ll)) remains form-invariant, 
S j W,j',\ , } = S j [i>,j,\} 



(35) 



For transformations with unit Jacobian this also holds 
for the partition function, Z[j', A'] = Z\j, A]. 

More generally, for arbitrary regular t the Jacobian of 
the functional measure is responsible for an anomaly 

Z[j', A'] = det(t)Z\j, A], W[f, A'] - W\j, A] + In det(t) 

(36) 

With G' a p = dW[f, \'}/df af3 related to G afj by 

G a p = t aa itpf}'G' a , p, (37) 

we find 

T B [G 1 , A'] = Tb [G, A] — In det (t) . (38) 
On the other hand, we note the transformation property 
det G = det(i 2 )dct G' (39) 

which implies 



1 



1 



T B [G',X'} --In det G' = F B [G,A]--m dot G (40) 

We may therefore write 

T B [G, A] = 1 In det G + f[G, A] (41) 

where the "reduced effective action" T is form-invariant 
under a simultaneous transformation of G and A accord- 
ing to eqs. (37) (34), i.e. f[G',A'] = f[G, A]. We con- 
clude that the BEA accounts for possible anomalies in a 
very simple and explicit form. 

The structural information on the possible form of Tb 
will be very useful for possible truncations. For example, 
let us exploit the invariance of T[G, A] under simultaneous 
transformations of G and A. This greatly restricts the 
possible independent invariants on which T can depend, 
i. e. 



G 



12 = —^^uPy8^a'f3'f'S'G aa 'GppiG 17 'Gs5' 

13 = -XafjyS^a'fl'y'S'GapGa'p'G^^iGsS' (42) 

These and other invariants involving higher powers of G 
are constructed by contracting each index of A with an 
index of G. We observe that 1% and 1\ are not 2PI. The 
three loop BEA must therefore be proportional to 12- 

Let us now come to the symmetries of T B . These are 
all transformations acting on ip which leave Tb invariant, 
with A kept fixed. Indeed, the full bosonic effective ac- 
tion rs[G] is invariant under all linear transformations 
(33) which obey det(t) = 1 and leave X a /3jS invariant, 
i. e. \ a 'i3'j'6' ta'a tfj'p ty y ts>s = A Q/ 3 7 5. This may 
be a very large symmetry group, much larger than the 
symmetry leaving Sj[ip] invariant for fixed A and j. For 

example, a pointlikc interaction A J d d x(ip(x)ip(x)) is 
invariant under local U(l) gauge rotations (with oppo- 
site phases for ip and ip) whereas the quadratic fermion 
kinetic term may not preserve the local symmetry. Nev- 
ertheless, the bosonic effective action Tg exhibits the full 
local symmetry. Actually, G(y,x) — (ip(y)ip(x)) trans- 
forms with different phases at x and y, i. e. G'(x,y) = 
exp i(a(x) — ot{y)) G(x,y), similar to a link variable or 
a string between x and y. We will give an example for 
the possible practical use of this large symmetry in ap- 
pendix A. For the pointlike interactions discussed in the 
next section the symmetry of Tb consists of local linear 
transformations with unit determinant, combined with 
appropriate space-transformations, e.g. the lattice sym- 
metries. 

Finally we mention the existence of useful transforma- 
tions which change the interaction. For example, mod- 
els with attractive and repulsive interactions are mapped 
onto each other by a transformation (33), (34) with the 
property X' a/3jS — —\ a p-y5- We present in sect. X a map 
between the repulsive and attractive Hubbard model. 
This establishes, for example, a direct link between an- 
tiferromagnctism in the repulsive Hubbard model and s- 
wave superconductivity and charge density waves in the 
attractive Hubbard model. 



V. LOCAL INTERACTIONS 

Several very interesting fcrmionic models assume a lo- 
cal four-fermion interaction. We concentrate here on a 
single spinor field (ip\(x), ip2(x)) = {4''[{ x )i V'j.( a; )) with 
an "antiparticle" (ip 3 (x) , ip4(x)) = (ipf(x),ip±(x)) (or 
electrons and holes). In a notation with 

*-(£)■*-(!) 

the action reads 



Sj = -\j a pMp + Jc{x) (44) 

with 

C{x)= l -\(^{x)iP{x)) 2 

= ^(^(^1(^+^2(^2(0;)) (45) 

In this section the integral J corresponds to a sum over 
points 1 in a D-dimensional lattice and includes a sum 
over discrete Euclidean time points r = em T , m T g Z, 
as appropriate for quantum statistics. The spinors are 
antiperiodic in the r-direction with periodicity given by 
the inverse temperature T _1 

M^r + T- 1 ) = -ip b {x,r) (46) 

We employ a,b — 1 ... 4 for the internal indices whereas 
a, [3 count all Grassmann variables, e. g. (3 = (6, x) — 
(6, r, x). The limit e — > has to be taken at the end 7 . 
More cxplicitely, one has 

* a = M*)= ( 47 ) 
and the local interaction corresponds to 

A Q( 3 7 5 = X a bcd(x,y,z,w) 

Xabcd^y,x^z,x^w,x > 

Xabcd = -Xe abcd (48) 

The interaction (45) has a high degree of symmetry 
which can be combined from independent local SL(4, C)- 
transformations among the four components ip a (x) at ev- 
ery point x and the symmetries of the lattice (i. e. appro- 
priate translations, rotations and reflections). The group 
SL(A, C) consists of all complex 4x4 matrices with unit 
determinant, i.e. t a p = t a b(x)S xy , det t a b = L There is 
no anomaly for this symmetry and the bosonic effective 
action is invariant under the corresponding transforma- 
tions of G a (3- Different models like the Hubbard model 
[1] or the Gross-Neveu [23] model can be obtained by 
choosing different sources j a p at the end. These sources 
will reduce the symmetry. Nevertheless, all these models 
will be described by the same bosonic effective action Tg ! 

Hermiticity of the Hamiltonion is reflected by 
Osterwalder-Schrader positivity [24] of the functional in- 
tegral. This means that the action Sj should be invariant 
under the transformation {-d 2 = I) 

z9(V(t,x)) = ^(-t,x) (49) 



7 In a continuum notation one has 
J x = e^a' 13 Jdrd D x,6 x , y = ea D 8{r - r')S D {x - y) with 
a the lattice distance. 
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if accompanied by complex conjugation of all coefficients 
and total reordering of the Grassmann variables, e.g. 

&(j>a(T,x)i> b (T',y)^ = 6 W 1pb'(-T' ,y)0 aa >-4> a >{-T,x) 

= -0 a a4a>(x')A>(y')6f, b (50) 

In a convenient basis with spin notation (47) one may 
have 8 



1 

1 



For "even sources" obeying 

jab{x,y) = -6l a ,j* a , b ,{x' ,y')6 vb 



(51) 



(52) 



the action Sj is indeed invariant and W[j] therefore real. 
Discarding spontaneous breaking of the ^-symmetry this 
implies 

#(G ab (x,y?) = -6 aa ,G a , v {x',y')6l b {x,y) 

= G* ab (x,y) (53) 

and we may restrict the bosonic fields G a p to those obey- 
ing eq. (53) without changing our previous constructions. 
For the choice (51), (53) one finds 

${x)$(x)) = $(x')i>(x'))* 
$(x)t^{x)) = ${x')t$(x'))* 
$ 1 (x)$ 2 (x)) = (54) 



Hubbard model 

In the following we concentrate on the Hubbard model 
while the Gross- Neveu model and its relation to the Hub- 
bard model are discussed in appendix A. The Hubbard 
model has a repulsive local interaction U > with action 
S given by 



(55) 



Here i, j denote the lattice sites in a D-dimensional cubic 
lattice. We concentrate on the simplest version where 
Tij — —t for next neighbours and zero otherwise. The 
identification with our previous discussion holds for A = 



In some other basis the transformation 8 can be appropri- 
ately modified. 



ell. A local source term is given by the chemical potential 
fj, which vanishes only for half filling 



ji 3 (x,x) = jn{x,x) = -fie 



(56) 



Another source term arises from the next neighbour in- 
teractions (or hopping term) 

j 13 (x , t;x±S,t) = j 2 i{x, r; x ± a, r) 
= -.?3i {x, t; x ± a, r) = - j 42 (x,t; x ± o, r) = -te (57) 

Here a is the unit lattice vector, a = \a\ the lattice dis- 
tance and x = (t,x),x = am with m a set of integers 
associated to i. Finally, the r-dcrivativc in J drtpd T ip is 
expressed by 

ji 3 (x , r; x, t + e) = -j\s{x, r; x , r - e) = 
J24,(x, t; x, t + e) = -]2a{x, t; x, t - e) = 
j 3 i (x, t; x, t + e) = -j3i (x, t;x,t-c) = 

342 (x, t; x, t + e) = -j i2 (x , r; x , r - e) = - - (58) 
One should take e — > at the end. 



VI. GAP EQUATION 

For local interactions the last term in eq. (23) is local. 
In the lowest order approximation the inverse fermionic 
propagator can therefore be written as the sum of the 
inverse classical propagator, — j, and a local "gap" A, 

G ab( x >y) = ~jab{x,y) + A ab (x)6 xy (59) 

where 5 xy is an appropriate generalization of the Kro- 
nccker symbol. Indeed, a nonzero A in certain chan- 
nels often induces an effective mass gap. In terms of the 
fermion propagator at equal arguments 



9ab(x) = G ab (x,x) 



(60) 



the generalized gap equation obtains by inversion of eq. 
(59) 



9ab{x) = {-j + A)J(x,x) 



(61) 



We further note the leading order relation (23) between 
A and g 



A ab (x) = -^e abcd g cd (x) 



(62) 



This turns eq. (61) into a closed "gap equation" for A 

Aab(x) - -^e abcd (-j + A)- d \x,x) (63) 

It is instructive to write this gap equation in Fourier 
space. We define 
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G ab (x,y)= [ [ e-^e^G ab (p,q) , 

J p J q 

G ab (p,q) = -G ba (-q,-p) (64) 
where the momentum integration reads 
d d p 



f _ D f d 1 



S pq = (ea D ) '(2n) d 5 d (p-q) 



(65) 



On a cubic lattice the integration interval is restricted 
by \pk\ < 7r/o and <5 D (p — 9) has to be taken modulo 
2ir/a. Finally, for r on a torus with circumference T -1 , 
as appropriate for the Matsubara formalism with p — 
2irnT, we take 



S pq = (ea D T) _1 S mn (2nf S D (p- q) 



d D p 



(66) 



The po-hitegration is bound by |po| < 7r/e or, equiva- 
lently, the Matsubara sum extends over a range \n\ < 
l/(2eT) with 5 mn taken modulo (eT)' 1 . (The largest 
possible value for e is T _1 .) The gap equation involves 
the momentum integration characteristic for a one loop 
expression 



9ab(Q) 



c g a b(x) 



= J (-3 + A)^( 9 -Q,q) = -9ba(Q) (67) 
We will concentrate on translation invariant sources 

jab(p,q) = Jab(q)S pq , J ab (q) = -Jba(-q) (68) 

A particularly simple situation arises for a translation 
invariant gap A ab (x) = A ab , A ab (p,q) = A ab 6 pq where 
(—j + A) is diagonal in momentum space and can there- 
fore easily be inverted. One finds that also g ab is trans- 
lation invariant, g ab {x) — g ab , and the homogeneous gap 
equation 



^ab = ~e abcd j (-J(q) + A)J 



(69) 



only needs the inversion of a 4x4 matrix for every q sep- 
arately. More generally, the gap can reflect spontaneous 
breaking of translation symmetry if A(Q) ^ for Q ^ 0, 
where Q = q — p and 



A(p,q) = A{q-p) = A(Q) 

e~ iQx A{x) = -A T (Q) 



(70) 



We will see in sect. VIII how to solve the gap equation 
(67) in the case of an inhomogeneous gap. 



VII. FREE ENERGY FOR LOCAL GAPS 

In general, the free energy functional (18) is a rather 
complicated object. We will show here that in the local 
gap approximation (corresponding to the leading order 
SD-equation) it can be reduced to a comparatively simple 
functional of the local gap A(x), namely 



Fi F 1 . 

Y = Y 2 Tr ln( ~- 7 



A) 



(71) 



where 



Fo 
T 



= ^=e abcd f A ab (x)A cd (x) (72) 
8A J x 



The first term i*o is a simple mass like term for the local 
field A, whereas the second term in eq. (71) reflects 
the fermionic fluctuation determinant in presence of the 
gap. The explicit expression (72) assumes the minimal 
local interaction discussed in sect. V. The construction 
given below can easily be generalized to more complex 
interactions. 

We recall at this point that a reconstruction of the 
free energy by integration of the gap equations (cf. [9]) 
is problematic. In general, one may have the knowl- 
edge about the gap equations in several directions Aj 
in the space of all possible A. This amounts to a 
set of partial differential equations which are equivalent 
to dF/dAi = 0. However, the integration of a given 
gap equation yields a function T(Aj) about which one 
only knows that it has a partial extremum in the Aj 
direction at a value Af 1 ^ which coincides with an ex- 
tremum of F. This is clearly insufficient for a compari- 
son of different local minima. The unsatisfactory situa- 
tion may be somewhat improved by combining gap equa- 
tions for more than one direction and by the inclusion of 
sources. Nevertheless, the available information is often 
even insufficient for a check of the intcgrability condition 
d 2 F/dA l dA J = d 2 F/dAjdAi. Our direct construction 
of the free energy functional avoids all these problems. 

Let us next derive the relations (71), (72) and discuss 
their applicability. For this purpose we treat the local gap 
approximation (59) as an ansatz for G _1 . In addition to 
the local field A we also use the local field g as defined by 
eq. (60). Since the lowest order relation (61) is not exact 
we treat A and g as independent variables. We will see 
below that the gap equation (61) results as an extremum 
condition for F in leading order. Our construction can 
be generalized beyond leading order where the relation 
between A and g becomes more complex. The free energy 



(73) 



obtains from eq. (25) and can be written in the form (71) 
with 



Ft 
T 



\ 1 f 

t = 2 Tr {Gj} + -Xabcd J gab(x)g c d(x) + c (74) 
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We next use j — — G 1 + A and the locality of A. Up to 
a shift in the irrelevant additive constant this yields 

F 1 f If 

Y = ~0 A °-b(x)gab(x) + ^Kbcd / 9 <ab{x) 9 cd(x) + C 

J X J X 

^ (75) 

The functional F\ depends on two independent local 
variables A and g as well as on j. We observe that the 
variation with respect to A at fixed g and j precisely 
yields the gap equation (61). (A similar procedure will 
fix A(g) also beyond the leading order.) Inserting A(g) 
the free energy F becomes a functional of g and one may 
then look for its minimum. For practical computations of 
the extrema it is more convenient to solve first the field 
equation for g as a functional of A (from the variation of 
Fi with respect to g at fixed A) and reinsert the solution 
g s [A] into eq. (75). From dF 1 /dg = one obtains 

-A ab (x) + ^Xabcdgcd(x) = (76) 

or, inserting (48) 

9ab(x) = --^e abcd A cd (x) (77) 

Reinserting (77) into the formula (75) for F yields the 
promised simple quadratic form (72). 

It can easily be checked that the variation of i*i (71) 
with respect to A(x) yields the local gap equation (63). 
The value of F\ at the minimum also coincides with the 
second formula in eq. (27). It is important, however, 
that the simple form (71) holds not only for the value 
of Fi at the minimum, but can also be used to search 
the extrema of Fi\ This turns the free energy formula 
(71) into a powerful tool for a comparison of different ex- 
trema of the free energy - as needed, for example, for the 
description of a first order phase transition. Neverthe- 
less, even in lowest order of the loop expansion for local 
gaps we have only shown the validity of the free energy 
(71), (72) for the location of the extrema and the values 
of F at the extrema. This is sufficient for a computa- 
tion of the ground state in lowest order, but not for the 
bosonic masses and succptibilities which involve second 
derivatives at the minimum. We believe, however, that 
our simple formula (71), (72) can serve as a good approx- 
imation for the bosonic masses and interactions in many 
practical situations. 

The free energy (71), (72) also establishes a close anal- 
ogy between the BEA and partial bosonization where 
the term quadratic in the bosonic field results from the 
Hubbard- Stratonovich transformation. In sect. XI we 
will discuss a definition of a free energy suitable for par- 
tial bosonization and make the comparison more explicit. 
In fact, the SD-approach and MFT lead to the same gap 
equation and the same free energy provided the partial 
bosonization is chosen precisely such that the quadratic 



term Fo coincides in both formulations. This provides for 
a very simple criterion for the possible choices of partial 
bosonization which are consistent with the Hartree-Fock 
approximation. 

We finally observe that the ansatz 

G a b{x,x) = g ab {x) , 

G ab (x,y^x) = (-j + A)- b 1 (x,y) (78) 

can be used for a simple estimate of the size of the higher 
order corrections in eq. (30) (31). The three loop contri- 
bution in eq. (31) has a local term ~ X 2 g 4 and a nonlocal 
term 9 <~ A 2 (— j + A)~ 4 . In particular, the size of the lo- 
cal correction can easily be compared with the leading 
order term ~ \g 2 once g has been computed in leading 
order. 



VIII. ANTIFERROMAGNETISM IN THE 
HUBBARD MODEL 

For the Hubbard model with repulsive coupling one 
does not expect a homogenous gap. For example, anti- 
ferromagnetic behavior corresponds to a flip of sign for 
neighboring lattice sites 

(ip(x)Tip(x)} = —(ip(x + a)Tip(x + a)) 

= —a exp(-iQax) (79) 
ha 

with Q a = (0, 7r/a, ir/a, . . . ). (We consider here a D- 
dimensional cubic lattice (d = D+l) with lattice distance 
a and have introduced the "Yukawa coupling" h a for later 
convenience.) For the condensate (79) the translations 
by two lattice distances 2a do not change the state. We 
therefore investigate gaps with the reduced symmetry of 
translations by 2a: 

A(Q) = A h S Qfi + A a S Q . Qa , 
A(p, q) = A h S p . q + A a 5 P+Qa , q . (80) 

Below the homogenous gap A^ will be connected to 
charge density and the inhomogeneous gap A a to an- 
tifcrromagnctism. 

With respect to the fundamental translations t x , t y by 
one lattice distance the gaps Ah — ~Aj t and A a = — Aj 
have even and odd parity whereas J is even 

t x (J) = J , t x {A h ) = A h , t x (A a ) - -A„. (81) 
This tells us immediately that for a similar decomposition 

(9h = -9h, 9a = ~9a) 

9(Q) = 9h $Q,o + 9a S Q , Qa (82) 



9 The part ~ (— j + A) 4 (x,x) has to be substracted in order 
to avoid double counting. 
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the "antiferromagnetic part" g a must involve an odd 
power of A a . For the Hubbard model with next neigh- 
bour interactions the nonvanishing sources read 



Ji3(-q) = J2i(-q) = -J3i(q) = -./42(a) 



D 



— [2irinT — p — 2t ^ cos(agfc)]e 



(83) 



k=i 



Our task is the inversion of the matrix (— j + A) in the 
gap equation (67) for the case of the inhomogeneous gap 
(80). 

Let us restrict the discussion to a gap 

(A h )i3 = (A h ) 24 = pe, (84) 
(A a )i 3 = -(A a ) 24 = -a 3 e, (A a ) 23 = -ode + ia 2 e, 
(A a )i4 = -cue - ia 2 e (85) 

where p transforms as a charge density and a = 
(0,1,0,2,0,3) is a real antiferromagnetic spin vector. We 
find that the pieces ~ a 1,0,3 commute with — J + Ah 
whereas the one ~ a 2 anticommutes. On the other hand, 
the contributions ~ 01,03 are real and the one ~ a 2 is 
purely imaginary. This implies for real hi , p 



A* a {-J{q) + A h f = {-J{q) + A h )A 



(86) 



The inversion of the matrix (— j + A) in cq. (67) can now 
be achieved by using the identity 

(-j + A)- 1 (p,q) = 

-M- 1 (p){ (j(p + Q a ) - A h ys pq + A* a 6 p+Qatq } (87) 
with 

M{ P ) = ( J(p + Q a ) - A fc ) * (j(p) - A,) - A: A a 

= ( 2 {A{p) + AirinT^s diag(l, 1, -1, -1)} 
A(p) = (27mT) 2 + a - fj, 2 eS + At 2 ( ^ cos ap^j (88) 



depending on 



Mcff = M + P i " 



(89) 



Here M 1 (p) plays the role of an effective squared 
fermion propagator 

M- 1 {p)=M- 1 {p)M*(p)e- i , 



M( P )=M*{p)M{p)e- A 

= (2nnT) 2 + a + tt 2 (J2 cos ap k ) 2 - fi 2 eS 
k 

+Apl s {2nnTf 



(90) 



For a > we observe that M (p) is strictly positive for 
all p, even for T = 0. An antiferromagnetic condensate 




FIG. 1. Phase boundary for the antiferromagnetic (AF) 
phase in the Hubbard model. We consider next neighbour 
interactions with U = t. The modifications of this leading 
order diagram due to bosonic fluctuations are important [13], 
see sect. XVI. 

therefore acts as mass gap for the fermions. Inserting our 
results in (67) one finds that g a b(Q) vanishes except for 

3(0) - - J^'Hq) (j(q + Qa) - a,,) * 
g(Qa) = I M-\ q )A* a 

Jq 

In leading order one has (A = Ue) 



g(Q a ) - U- 1 



a k T k 
-a k Tl 



This yields the antiferromagnetic gap equation 
U 



■ f M-\q)A{q) = 1 

J q 



(92) 



(93) 



where we recall J q = ea D TJ2(^)- D f d D q, q\ < (ir/a) 2 

n 

and n is half integer. As it should be, this gap equation 
becomes independent of e 



a D UTj2 



d D q 



A(q) 



(2tt) d A 2 (q) + (4irnTp eS f 



(94) 



The solution of the gap equation (94) determines the an- 
tiferromagnetic order parameter a as a function of T and 
p c ff. In particular, it allows us to determine the region 
in T and p, for which the antiferromagnetic state {a > 0) 
occurs. We show this region in fig. 1. By construction 
the gap equation (94) is one loop exact and therefore 
consistent with Hartrcc-Fock [25]. Nevertheless, we have 
implemented the spin rotation symmetry in a simple and 
direct way and can therefore proceed to a straightforward 
computation of the free energy for the antiferromagnetic 
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state (below) and the spin wave excitations above the 
antifcrromagnetic ground state [26,13]. 

Eq. (94) coincides with the gap equation in the 
mean field approximation derived in the "colored Hub- 
bard model" [7] if we identify (cf. eq. (79)) a = 
h 2 a — h 2 aa, p, e f[ = h p p, U = 2h 2 /n 2 . In con- 
trast to the mean field approximation where the rela- 
tion between h 2 and U depends on free parameters Aj 
which reflect the details of the bosonization procedure 
(hi = (2Ai + A 2 -3A 3 + l)7r 2 J7/3), the present Schwinger- 
Dyson equation gives a unique relation between the four 
fcrmion coupling U and the effective Yukawa coupling 
h a . (The latter characterizes the interaction between the 
composite boson fields a(x) responsible for antiferromag- 
netism and the fermions.) In particular, the part of the 
phase diagram shown in fig. 1 corresponds to one of 
the phase diagrams shown in [7], now for a particular 
value h\ = (ir 2 /2)U. For half filling (p = 0),t = U and 
U « leV the transition temperature is well compatible 
with the observed range of transition temperatures for 
the antifcrromagnetic phase in materials which exhibit 
high temperature superconductivity in the case of doping 
(e.g. /j ^ 0). The inclusion of the bosonic fluctuations 
will lower the critical temperature (sect. XVI). 

For the free energy of the antifcrromagnetic state we 
can directly use the formula (72) and insert eqs.(84)(85), 
resulting in 

Y = hS dT ^ {a -~ p2) (95) 

i 

The term ~ a coincides with the result of partial 
bosonization [7] precisely for the choice h 2 = (ir 2 /2)U. 
For given p e g the computation of the free energy can 
therefore be directly inferred from [7]. On the other hand, 
the negative sign of the term ~ p 2 indicates that there 
exists no consistent choice of partial bosonization which 
reproduces the gap equation of BEA for real p. (A consis- 
tent choice must have a positive quadratic term.) This 
demonstrates in a simple fashion that it is not always 
possible to define a well defined Hubbard-Stratonovich 
transformation such that the mean field approximation 
is consistent with the Hartrec-Fock result! 

A similar investigation can be done for the possibility 
of superconductivity. A local d,.2_ y 2-wave superconduct- 
ing order parameter reads 

TT 2 2 

i(ip(x)T 2 ip(x)) = 7rr d Yl exp(-iQ^fc) 
ztld k=i 

TT 2 - 2 

i{iP(x)T 2 <P(x)) = -— dJ2 expi-iQ&k) (96) 
with Q\ = (0, tt/o, 0, 0, . . . ), Q 2 = (0, 0, n/a, 0, . . . ). The 



discrete symmetry ■& (cf. eqs (49), (51) implies 10 d = d* . 
Similarly, a d^-wave or s-wave superconducting order 
parameter takes the form 

2 

i(ip(x)T2ip(x)) = ^-c exp ( - i(Q a x)j (97) 

or 

TT 2 

i(^(x)T2ip(x)) = —s (98) 

llg 

Here IT2 assures the antisymmetry and all above order 
parameters are invariant under spin rotations. 

In particular, for a pure d x i _ y 2 -superconductor only d 
differs from zero and we may consider a local gap 

I 2 

A(p,q) = A h S p , q + -A d ^d p+QU (99) 
fe=i 

with A h given by eq.(84) and (A d )i 2 = d, (A d ) 34 = -d* 
or 

The computation can now be performed in analogy to 
the antiferromagnetic gap. For the local gap (100) the 
"classical free energy" reads 

Y = ~2ljj dT ^* (101) 

i 

Due to the negative sign in eq. (101) one finds no solution 
for the lowest order gap equation with d*d ^ 0. 

Nonlocal superconducting gaps are well motivated by 
functional renormalization group studies [20-22]. They 
have been investigated in the framework of partial 
bosonization in [7]. Unfortunately, the Fierz ambiguity 
influences strongly the MFT phase diagram. A compu- 
tation of the BEA would therefore be quite useful. How- 
ever, in lowest order all nonlocal superconducting gaps 
vanish, i.e. G\2(x,y ^ x) = 0, in the case of local in- 
teractions. A nontrivial gap equation for a nonlocal su- 
perconducting order parameter must therefore proceed 
to the next to leading order equation based on eqs. (22), 
(30). As an alternative, we notice that the renormaliza- 
tion flow of the fermionic interactions induces nonlocal 
interactions once the fluctuations with q 2 > k 2 are in- 
cluded. This is precisely the aim of the renormalized gap 
equation derived in sect. XIII. Within the BEA the inte- 
gration of the high momentum fluctuations with q 2 > k 2 
can be performed by an approximate solution of the exact 
functional renormalization group equations which will be 
derived in sect. XII. A proposal for the stabilization of 
d-wave super conductivity in an approach similar to the 
BEA can be found in [27]. 



This could be different for other realizations of 1?. Replac- 
ing the unit in eq. (51) by T3 would lead to d = —d*. 
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IX. SYMMETRIES OF THE HUBBARD MODEL 

In this section we consider more closely the symmetries 
of the Hubbard model. They will give important con- 
straints for the general form of the free energy. We will 
restrict the discussion here to those transformations that 
also leave the nonlocal source terms (57) (58) invariant. 
On the other hand, the chemical potential (local source 
(56)) will only appear through the field equation as a 
possible symmetry breaking term. Besides the discrete 
lattice symmetries the next- neighbour Hubbard model 
exhibits a global 5(7(2) x SU{2) symmetry. 

Consider first local infinitesimal transformations 



Sip a (x) = iL ab (x)ip b (x) 



(102) 



The interaction term is invariant provided TrL = and 
invariance of the term in eq. (55) involving d T , 



requires 



L K + KL = 



(103) 



(104) 



Both conditions are obeyed of L generates the symmetry 
517(2) x 5(7(2), with 



L{x) = a(x)S + 0(x)U 



where 5 are the spin generators 



Si = 







(105) 



(106) 



and 



r 2 
t 2 

1 

-1 



U '=2 





«T 2 



-IT2 





(107) 



obey the commutation relation [Ui,Uj] = ieijkUk, 
[Si,Uj] = 0. The electric charge generator is given by 
Q = —2Us whereas U\ and (7 2 change electrons into holes. 
We note that the chemical potential respects the spin and 
charge symmetry 5(7(2) x (7(1) while it is not invariant 
under the transformations generated by U\ and (7 2 . 

For the Hubbard model on a quadratic (cubic) lattice 
with next neighbour interactions we may divide the lat- 
tice into two sublattices A, B, such that the hopping term 
~ t connects spinors on different sublattices ~ ipA^B- 
The next neighbour interaction remains invariant under 
a global 5(7(2) x 5(7(2) symmetry with 



OLA = as, 03A = 03B,P 



l A 



-01B, 02A = -02B (108) 



The Hubbard model with next neighbour interaction has 
no further continuous symmetry beyond 5(7(2) x 5(7(2). 



Adding diagonal interactions reduces this symmetry to 
5(7(2) x (7(1). 

The spinors t/j a transform as the (2, 2) representation 
of 5(7(2) x 5(7(2), which takes its canonical form in the 
basis 



tpi , -itpi 

-■03 , -#2 



Equivalently, we can also use 



= I(af0-/30f) (109) 



ltp 3 , I1p4 ' 



(110) 

In consequence, the spinor bilinears g a b belong to two 
irreducible representations (3,1) + (1,3). The magneti- 
zation 

m(x) = (0(x)t0(x)) = (i> a (x)(KS) ab Mx)) ("1) 

is a vector with respect to the spin rotations and a singlet 
with respect to the generators U. On the other hand, 



u(x) = (^ a (x)(KU)abipb(x)) 



(112) 



is a vector with respect to the (7-rotations and a spin 
singlet. In particular, an antiferromagnctic state corre- 
sponds to m(x) ~ acxp(— iQ a x) and the charge density 
p{x) is given by p(x) = Uz{x) — (ip(x)?p(x)). The six local 
fcrmion bilinears g a b(x) can be written as linear combina- 
tions of rrii{x) and Ui(x). In particular, the electron-pairs 
(or hole pairs) 

tt(x) = i(ip(x)T2ip(x)) = 2ip 1 (x)ip2(x) , 
ir*(x) = -i(^(x)T 2 ^(x)) = -2^ 3 (x)Mx) (113) 

can be written as 



7T = U 2 + IU\ , 7T = M 2 — IU\ 



(114) 



Consistently with cq. (53) the fields m and u are real. As 
a consequence of the 5(7(2) x 5(7 (2) symmetry the effec- 
tive action can only depend on invariants as Uk(x)uk(x) 
etc. In particular, the free energy (95) can be extended 
to 



Y = Tij dT ^ a ~ ^ ' 7 = ~ p2 + 7r * 7r 



UfeUfe (115) 



(for homogeneous tt ~ s). For \i ^ 0, the degeneracy 
between the three components of u will be lifted by the 
solution to the field equation, since the source ~ [i points 
in the three-direction and therefore breaks the invariance 
with respect to [/-transformations. 



X. MAPPING BETWEEN REPULSIVE AND 
ATTRACTIVE HUBBARD MODEL 

In this section we discuss a mapping from the repulsive 
(U > 0) to the attractive (U < 0) Hubbard model. The 
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automorphism 

T A :<f)->i4> , -> -i0 (116) 

maps the generators S and U onto each other. (On the 

level of individual components one has 

T A : 4 i0 2 , 02 -> ~«04 , ipi->ipi , ip3 -> ip3-) 

It changes the sign of the interaction term (?7 — > — [/) and 

leaves the term C T (103) invariant. We may consider a 

second automorphism 

T B : — > -iT 3 4>T 3 , 4> — » «T 3 0T 3 (117) 

that also changes the sign of the interaction and leaves 
£ T invariant. (It reads in components Tg : 04 — > 
— i02, 02 - * £04, 01 - * 01, 03 — > V"3-) This transfor- 
mation differs from Ta by an additional flip of sign of 
the "lower spin components" 02, 04- The next neighbour 
coupling of the Hubbard model remains invariant if T4 
acts on the sublattice i\>a and Tb acts on 0s . For p = 
we have therefore established a transformation that maps 
the attractive on the repulsive Hubbard model and vice 
versa. 

In order to investigate the interesting consequences of 
this map we should have a look at the action on the local 
fermion bilinears m and u. For the sites of the sublattices 
A and B one has different transformation properties 

rh — > u , u — > rh for A 

rh — > u , u — > m fori? (H8) 

with ui ; 2 = — Ui,2)U3 = "3 and similar for rh. In conse- 
quence, the antiferromagnetic state of the repulsive Hub- 
bard model is mapped into a charge density wave with 
(0(x)0(x)) ~ (— l) p |a| in the attractive Hubbard model. 
(Here P = 1 on A and P = — 1 on B.) It is convenient 
to introduce the field v(x) as 

v 3 (x) = u 3 (x) , vi t2 (x) = exp(iQ a x)ui t2 (x), (119) 

which obeys 

. f u(x) for A , 

= { 4) for 5 ( 12 °) 

In terms of this field the combined automorphism T (i.e. 
f = T A on A, f = T s on B) simply maps 

f : m(a;) ^ w(a;) (121) 

With respect to the global {/-transformations (108) the 
field v(x) indeed transforms as a global vector (i.e. the 
same S'[/(2)-transformation for every x), in complete 
analogy to m(x) with respect to the 5-transformations. 
(For the transformations U\,U2 the minus sign in eq. 
(108) for B is absorbed by the minus in the definition of 

«1,2-) 

The mapping T (121) allows a simple interpretation 
of several properties of the Hubbard model. For half 



filling (ji — 0) the low temperature state of the repulsive 
Hubbard model is antiferromagnetic. Without loss of 
generality we may choose the 3-direction (m(x)} = m 3 (x) 

m 3 (x) = a 3 enp(iQ a x). (122) 

This implies, for the same low temperature T and the 
same \U\, that the attractive Hubbard model is charac- 
terized by a charge density wave 

p{x) = pexp(iQ a x) , p = a 3 (123) 

The S'?7(2)-spin-symmetry of the repulsive Hubbard 
model exhibits a degeneracy of the order parameter in 
the (ai, 02, (i3)-space. Correspondingly, the [/-symmetry 
of the attractive Hubbard model predicts a degeneracy in 
the space (vi,v 2 ,v 3 ). The 5'C/(2)-partners of the charge 
density wave (123) are the homogeneous s-wave Cooper- 
pairs 

tt(x) = 7f , 7f*7f = aj (124) 

Thus the low temperature state can also be s-wave 
super-conductivity. This degeneracy between the charge- 
density wave and the s-wave superconducting state for 
p, = seems compatible with existing numerical results 
[28], [29]. Furthermore, we predict for p, = the ex- 
istence of two (real) massless Goldstone bosons for the 
ordered state, both for U > and U < 0. 

A recent renormalization group study [13] has com- 
puted the temperature dependence of the antiferromag- 
netic order parameter a for the repulsive Hubbard model 
in two dimensions. These results can be taken over one 
to one for the temperature dependence of p or \n\ for the 
attractive Hubbard model. The mapping T can also be 
applied for the discussion of the correlation length £ for 
T > T c . Close to the critical temperature T c one finds 
[12] 

£«cexp(Wp) (125) 

Here T c depends on the size of the experimental probe 
and eq. (125) holds for a size of 1 cm. The compu- 
tation [13] yields for \U\/t = 3 a value T c w 0.115*. 
However, the value of T c is not universal and depends 
on details of the model (and the computation). Con- 
sistent with the Mermin- Wagner theorem T c vanishes 
logarithmically in the limit where the size of the probe 
becomes infinite. The same behaviour of the correla- 
tion length governs the (anti)fcrromagnetic correlation 
function (m(x)m(y)) exp i (Q a {y — x)) for U > and 
the charge density wave and superconducting correlation 
functions (p(x)p(y)) expi (Q a (y - x)) or (it* (x)ir(y)) for 
U < 0. 

For the attractive Hubbard model away from half fill- 
ing (p y£ 0) one expects that the degeneracy is lifted 
and the low temperature state is s-wave superconduct- 
ing, with (ir(x)) = 7T 7^ 0. The phase transition for d = 2 
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is expected to be of the Kosterlitz-Thouless [30] type, 
such that in the limit of an infinite sample only a suit- 
able renormalizcd order parameter differs from zero for 
T < T c [31]. Indications for this transition have been 
seen in simulations [28], [29]. These properties can be 
mapped to the repulsive Hubbard model in a homoge- 
neous magnetic field. In fact, the transformation T maps 
a source term with a chemical potential into one with a 
constant magnetic field: 

t : fj,p(x) = fj,v 3 (x) <-> B 3 m 3 (x) , B 3 = fj, (126) 

We conclude that the low temperature state of the repul- 
sive Hubbard model in an homogeneous magnetic field B 
(for fj, = 0) is antifcrromagnetic with direction orthogonal 
to B, i.e. for B = (0, 0, B 3 ) one has a 1>2 ^ 0, a 3 = 0. The 
spin symmetry is now reduced to the [/(l)-spin-rotations 
around the three axis (which is mapped by T to the 
charge-symmetry for U < 0). In complete analogy to the 
superfluid for U < one predicts for U > and B 3 ^ 
the characteristic properties of a spontaneously broken 
[/(l)-symmetry, namely the existence of one Goldstone 
boson and of stable vortex solutions. 

The T-transformation offers many interesting perspec- 
tives. Results from numerical simulations of the attrac- 
tive Hubbard model can be directly taken over to ap- 
propriate quantities in the repulsive Hubbard model. In 
the other direction, the presumed d-wave superconduc- 
tivity for the repulsive Hubbard model sufficiently away 
from half filling can be mapped into a corresponding spin 
wave in the attractive Hubbard model in an homogeneous 
magnetic field. 



XI. PARTIAL BOSONIZATION 

The BEA is economical in the sense that redundancy 
is avoided. In principle, Tb contains the complete infor- 
mation about all Greens functions. Many physical situa- 
tions, however, are well approximated by the coexistence 
of effective local bosonic and fermionic degrees of free- 
dom. These aspects of locality are not very apparent 
in the BEA, the latter being based on the bilocal fields 
G a p. For situations dominated by local bosonic degrees 
of freedom - as often for the critical behavior near a sec- 
ond order phase transition - an effective locality can be 
recovered by restricting the BEA to local fields g a b{x). In 
contrast, if local fermions are relevant, it may be advan- 
tageous to keep explicitely the local fermion fields ip{x) 
as well. 



Mixed effective action 

We may construct a "mixed effective action" Tm which 
depends both on ip a and G a f), thereby accepting a cer- 
tain redundancy of the description. This task is eas- 
ily achieved by a Legendre transformation of W^j] (cf. 



sect. II) with respect to both the fermionic and bosonic 
sources 

T M [i>, G] = -W[r], j] +Va^a + \j a pG ati (127) 



with 



^ = ^M,G^,i] = ^M (128) 



dr) a 

The field equations read now 



dj, 



a/3 



dT M _ ,. d£M_ _ . 

o / Va i — Ja/3 



(129) 



and the relation W^T^ = 1 is now a matrix relation 



in the combined space of fermionic and bosonic fields. 
The BEA is simply related to Tm by restriction to 

V = o, 



r B [G] =r M [V = o,G] 



(130) 



This follows from the fermionic character of ip which im- 
plies that Tm can only involve even powers of ip. In turn, 
T] a = —dTM/dipa vanishes for ip a = and the definition 
(127) coincides for with eq. (12). Furthermore, the mixed 
derivatives d 2 T m / '(dipdG) vanish for ip — 0. The matrix 
of second functional derivatives T M ' therefore becomes 
block diagonal in fermion and boson space if ip = 0, and 
similar for if rj = 0. We emphasize that the sim- 

ple relation (130) is closely connected to the fermionic 
character of tp and bosonic character of G. It would not 
hold automatically for the two-particle irreducible effec- 
tive action for a fundamental bosonic theory. We also 
point out that the fermionic effective action rjr[^] does 
not obtain for G = 0, but rather obeys 

r F [^;j] = r M [^G [V;j]] -j-G°[t,j} (i3i) 

Thereby G [ip;j] is a functional of ip which is determined 
by the solution of the field equation dTM/dG = j in 
presence of arbitrary tjj. With 

dT F 8T M 



\G 



one infers 



8 2 T F d 2 T M d 2 T M 3G°JV] 



dip a dipp dip a dip dtppdG-fS dip Q 



(132) 



(133) 



(For ip = the second term on the r.h.s of eq. (133) 
vanishes.) Taking two more derivatives and evaluating 
at ip = relates the terms quartic in ip in Tp and Tm 



d 4 T, 



M 



dip e dtp n dtp a dipp dip e dip n dip a dtp(3 

| d 3 T M d 2 G° s m _ d 3 T M d 2 G°M 
dxpadippdG^g dtpedipr! diprjdippdGjg dip e dip a 
d 3 T M d'G ( >M 



+ 



dipedippdGjs dipedipa 



(134) 
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The last three terms on the r.h.s. of eq. (134) corre- 
spond to "trees" from the exchange of G-bosons, as can 
be verified by taking suitable derivatives of eq. (128). 

The redundancy of the mixed description is reflected 
by exact identities for the derivatives of Tm- The relation 
(4) holds for arbitrary j and r\. It results in a functional 
identity for Tm which holds for arbitrary ip and G 



\ t m)ff) = G af} - tpaipp 

/ afj 



(135) 



-l 



Here (t m ^ denotes the fermion-fcrmion component 

of the inverse of the matrix of second derivatives T M '. 
For ip — one has (cf. eq. (133)) 

(r$)^ [V> = 0, G] = (if) 1 [V = 0; j] (136) 

and we recover from eq. (132) the relation (15). Further 
identities can be obtained by taking partial derivatives 
of eq. (135) with respect to ip or G. The identity 



M 



dip a dipp 



\i>=o 



a/3 



(137) 



can be used for a direct evaluation of G -1 without the 
need to solve a field equation. 

The redundancy of Tm can be greatly reduced by 
resticting the most general sources j a /3 to a small sub- 
set. For example, we may only consider local sources 
jab(x) and the associated local bosonic fields g a b(x) = 
SW/Sjab(x). Performing the Legendre transform similar 
to eq. (127) one defines an effective action Tfb which 
only depends on local fermionic and bosonic fields 



(138) 



T FB = -W + j ( T] a 1p a + -iab9ab 



with field equation 



dThU -Va{x) , ^-^-=j ab (x) (139) 



8ip a (x) 



Sg a b(x) 



The relation between fps and Tm proceeds by split- 
ting 



G a fi = 9ab{x)5 xy + G ab (x, y) 



(140) 



Next one solves the field equations for G ab (x, y) as a func- 
tional of arbitrary local fields g ab {x) 1 ip a (x). Inserting this 
solution into Tm yields Tfb- Using the redundancy of 
Tm we may also compute G ab (x,y) from the fermionic 
part of Tm 

G ab (x,y)= ((t^Y 1 ) (x,y)-g ab S xy (141) 

V V /FF /ab 

As an example, we can compute the two loop effective 
action Tfb from eq. (31). Using a local interaction one 
has 



f fb,2 = \ j Kbcd9ab(x)g c d(x) - ^Trln (Tm) Fp 



1 

~48 



^■abed^efgh 



(r£>)^) c ^((T^Jx^) (142) 



(2) 

with T M ' depending on g and We note that the full 
information about all n-point functions is contained in 
Tfb due to its dependence on the fermionic fields. The 
restriction to the bosonic part by setting ip = is now a 
functional of the local bosonic fields g a b(x) only. 



Local fermion-boson-model 

The partially bosonized effective action Tp B can be 
equivalently computed from a functional integral involv- 
ing fermionic and bosonic degrees of freedom. This be- 
comes apparent by a generalized Hubbard-Stratonovich 
[5] transformation. We restrict the discussion here to a 
local four fermion interaction (47) and rewrite the parti- 
tion function in an equivalent form as 



bed 



Z[V, J] = J exp j - jT ( A £a 

(j}ab{x) ~ ^ a (x)A{x) - jab(x)j 

(g C d(x) - i> c {x)ipd(x) - jcd(x)j 



-Is 

2^ °abcd 



-Va(x)^ a (x) 



(143) 



Here the totally antisymmetric part of S a bcd vanishes - 
otherwise S a bcd = S c dab = -Sbacd = -S ab dc is arbitrary, 
reflecting the "Ficrz ambiguity" [8,18]. The integral over 
g is a well defined Gaussian integral. It only yields an 
irrelevant multiplicative constant provided that 



£ m = --^Kbcdgab(x)g C d(x), 



^abed — ~ ^ e a 



bed 



s a 



bed 



(144) 



is a positive definite quadratic form. This can be estab- 
lished by the use of appropriate involutions acting on the 
■0 as complex conjugation. This involution relates the 
real and imaginary parts of g and correspondingly of g. 

The "insertion of a unit" has been chosen such that the 
terms quartic in ip cancel the quartic interaction in Sj [ip] . 
It is obvious that this cancellation does not depend on the 
choice of S a bcd - different choices of S a b c d "distribute" the 
original four fermion interaction into different channels of 
exchange of local bosons. The remaining fermionic action 
is quadratic in ip 
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/ 



Z = / ThpVg cxp {-(S m + S y + S F . km + Sj + S v )} 



1 r 



Sy = —Xabcd J 1pa(x)->Pb(x)gcd(x), 
S F ,kin = \ J i>a(x)'4>b(x)J a b(x) 



3ab{x,y)^a{x)-lpb{y) 



ijj a (x)K ab (x,y)ip b (y), 



- / T1 a (x)lp a (x) 

J X 



Here we have introduced the sources 



Jab(x) = --Kbcdjcd(x) 



(145) 



(146) 



which appear in appropriate linear combinations in eq. 
(143) 

jab(x) = PabcdJcd(x) (147) 

The "bosonic" pieces of the action are then given by 



Sj = - J (^Jab(x)g a b{x) ~ ^PabcdJab{x)Jcd(x) S j (148) 

According to our condition on A the matrix p a b.cd is pos- 
itive definite such that after diagonalization all "mass 
eigenvalues" in C m are positive. 

For a suitable choice of the local sources J a b{x) one 
may cancel the local part of the bilinear source j a b{x, y), 
e.g. mass terms for the fermions or a chemical poten- 
tial. These pieces are now reflected by terms linear in g. 
The nonlocal part of j a b(x,y), i.e. the kinetic operator, 
remains, however, and we end with a kinetic operator 



K ab (x,y) = -jab(x,y) + J a b(x)5 xy 



(149) 



We recognize in eq. (145) the standard formulation of 
a partition function with bosonic and fermionic fields, 
coupled by a Yukawa interaction S y . We emphasize 
that both the Yukawa couplings in S y and the bosonic 
quadratic term ( "boson mass term" ) in S m depend on the 
choice of S a bcd- We know from our construction that all 
exact results for correlation functions must be indepen- 
dent of Sabcd- On the other hand, simple approximation 
schemes like mean field theory (which integrates out the 
fcrmion fluctuations for a given local "background field" 
g) typically show a strong dependence on S a bcd- This 
often reflects a severe shortcoming dubbed "Fierz ambi- 
guity" [8] . The validity of an approximation can be (par- 
tially) tested by establishing the independence of results 

On Sabcd- 



In complete analogy with the standard treatment we 
may define a free energy Wfb as a functional of the local 
sources J a b{x) and r] a (x), while keeping K a b(x,y) as a 
fixed part of the action 

W FB [v, J] = ^Z[ri, J] + ^p a bcd J J a b{x)J c d{x) (150) 

Defining the effective action Tfb by the usual Legendre 
transformation 

IVb[V>,<7] = -W FB + J (Va^a + ^Jabfjab), 

_ SWfb , , N 5Wfb 

9ab{X) = 



SWfb . , , 

OJab(X) 



Sr] a {x) 



(151) 



we recognize that the only difference between Tfb and 
Tfb (cf- eq. (138)) results from the term quadratic in J 
in eq. (146) and the difference between J a b§ab and j a b9ab- 
The precise definition of j a b in eq. (138) and J ab in 
eqs. (145), (146) depends on the splitting of j a b{x,y) 
into a part kept fixed (K ab (x, y)) and a part with respect 
to which the Legendre transform is performed. Indeed, 
the redundancy in the construction of the effective action 
Tfb allows us to put a local fermion-bilincar cither into 
the source for g ab or into a part of the action. Let us in- 
corporate the kinetic term K ab (x, y) (149) into the action 
such that the source term in the second equation (139) 
becomes j a b(x) — > J a b(x). Comparing d\nZ/dj ab = g ab 
and eqs. (146), (148) yields then the simple relation 



lab{x) = (g a b(x)) = 9ab{x) + ^PabcdJcd(x) 
= {tpa(x)tpb(x)) + ^PabcdJcd(x) 



(152) 



This relates the fermion-boson model introduced via the 
generalized Hubbard-Stratonovich transformation to the 
partial bosonization from the BEA 

1 



Tfb = Tfb + -^Pabcd J J a b{x)J c d{x) 

1 f STfb STfb 

= T F b + ^Pabcd I t 7 r ~? -r-: (153) 

8 Jx ogab{x) 6g c d(x) 

Both Tfb and Tfb obey the same field equation 



FB 



ST 



FB 



5g ab (x) Sg a b(x) 



Jab(x) 



(154) 



If we interprete g a b(x) as an extremum of Tj t pB = 
Tfb - \ J x Jabfjab (at fixed J) then g a b is a corresponding 

extremum of Tj fb = T F b - \ J x J a b9ab- 

Away from the extrema, however, the functionals 
Tj^fb and Tjfb differ. For example, the second func- 
tional derivatives at the extrema do not coincide. This 
has important consequences for the mass matrices of the 
composite bosons g or g. Indeed, one obtains 
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FB 



S 2 T 



FB 



Sg a b(x)Sg cd (y) 



zPefgh 



I 



Sg a b(x)Sg cd (y) 

d 2 T FB d 2 f FB 
dg ab (x)dg ef (z) dg gh (z)dg cd (y) 



(155) 



by differentiating (dT FB /dg)(g + pJ) = (dT FB /dg)(g) = 
J(g),eq. (154), with respect to g. The difference between 
the two point functions for the bosons g and g simply 
reflects the fact that the choice of bosonic fields is not 
unique. The bosonic field g is a nonlinear function of the 
fermion bilinear g according to eq. (152), i.e. 



g a b{x) = g a b{x) 



1 5T FB 

Pabcd 



■led 



(x) 



[.9] 



(156) 



Of course, the computation of the correlation functions 
for g and g yields the same results if we either start from 
the mixed effective action or from partial bosonization 
via a Hubbard-Stratonovich transformation. 

We conclude this section with a simple but impor- 
tant observation. In general, the functional r^sls 1 ] is 
not bounded from below. On the other hand, the func- 
tional IVb[<7] contains a positive definite "classical part" 
inherited from the Hubbard-Stratonovich transformation 
(cf. eqs. (144), (145)) 



r (cl) _ 1 t 



bed 



g a b{x)g c d{x) 



(157) 



In many cases the fluctuation effects become unimpor- 
tant for large values of \g\. Then T F b grows quadrati- 
cally for large \g\. This is reflected by the positive def- 
inite second term in eq. (153). We also note that for 
any logarithmic divergence Y F b ~ In g, where T F b may 
get arbitrarily large negative values, the positive second 
term in eq.(153) will dominate such that T FB gets large 
positive values. We conclude that in contrast to T FB the 
functionals T F b arL d Tj iF b are typically bounded from 
below. The boundedness of Tj jF b may be helpful for the 
classification of extrema. 

In particular, we may identify Tj^ FB as a suitable defi- 
nition of a free energy. The solutions of the field equation 
correspond to its extrema. They are stable if the matrix 
of second derivatives is positive semidefinite. As we have 
argued Tj iF b is typically bounded from below. On the 
other hand its precise form depends on the particular 
choice of S ab cd (144). This corresponds to the depen- 
dence of the variable g on S a bcd according to eq. (156). 
As mentioned already before, the independence of the 
physical quantities on the choice of S ab cd can be used as 
a test of the validity of approximations. In sect. VII we 
have formulated a simple criteria for the choice of S a bcd 
which ensures that the solution of the lowest order gap 
equation agrees with the Hatree-Fock-approximation. As 
shown explicitely in sect. VIII this can be achieved at 
best for some of the directions in field space. 



XII. EXACT FERMIONIC RENORMALIZATION 
GROUP EQUATION 

So far we have used the Schwingcr-Dyson equation as 
a main tool for a computation of the bosonic effective 
action. This equation does not reduce to a closed set of 
equations when we go beyond the leading order gap equa- 
tion. The term ~ Y a p in eq. (22) involves the bosonic 

propagator (r^) -1 and we need, in turn, an equation for 
this propagator. It is obvious that the correction term 
~ Y a p becomes important whenever a bosonic bound 

state becomes (almost) massless. In this case T B ' has 

a zero eigenvalue and the bosonic propagator (r^) -1 
becomes very large in some range of momenta. In partic- 
ular, the lowest order Schwinger-Dyson equation is not 
expected to give a reliable description of critical phenom- 
ena associated to a second order phase transition. 

Furthermore, the expansion of Tb in powers of the 
coupling involves the microscopic or "bare" coupling A. 
This is not very satisfying from a renormalization group 
viewpoint. In fact, it should be possible to integrate out a 
certain range of short distance fluctuations k 2 < q 2 < A 2 
such that a new action Sj takes into account the effects 
of these fluctuations. As a result, an effective coupling A& 
will replace A. One therefore would like to derive some 
type of "renormalized gap equation" which involves Afe 
rather than A. In this and the next sections we will show 
how both problems can be attacked by the formulation 
of an exact renormalization group equation for the BEA. 

We will first base our investigation on the successful 
exact flow equation for the fermionic effective average 
action and translate it from the original formulation in 
terms of fermion fields to the bosonic language used in 
this work. For this purpose we introduce a fermionic cut- 
off such that only modes with momenta q 2 > k 2 are 
effectively included in the functional integral. (The mo- 
mentum q may measure the distance from the Fermi sur- 
face.) Our formalism will result in an exact renormalized 
gap equation which takes precisely the form of the lowest 
order equation (23). Only the couplings and sources are 
replaced by renormalized couplings and sources A& and 
3k 



T af3 — -\3k)af3 + ^(^afJjsG^S 



For k — > the cutoff term vanishes. 

We will also derive from the exact renormalization 
group equation for the BEA an approximate formulae 
for the ^-dependence of the renormalized coupling A& 

{(Afe)a 7 e J u(Afe) / g^ I/ - (Afc) < a 5en(^k) /3jr)i>} 

(159) 



R 



a/3 



(158) 



Gy,vG ep (dkR k ') G 

\ /pa 



as well as for the renormalized source term 
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G jS G^G ep d k (r[ f) ) G av (160) 

V / pa 

On the r.h.s. the vertices Xk are contracted with fermion 
propagators G evaluated at the minimum of Tj, i.e. so- 
lutions of cq. (158). The first equation therefore takes a 
one loop and the second a two loop (setting sun diagram) 

(F} 

form. The fc-dependence arises through the cutoff R k . 
The system of equations (158-160) is closed and can be 
used as a practical tool. A more explicite contact to criti- 
cal phenomena will be established in sect. XV by varying 
an infrared cutoff for the composite bosonic fluctuations. 

In this section we make use of the functional rcnormal- 
ization group for the effective average action for fermionic 
systems [10]. We recapitulate the main features here 
briefly. In order to derive an exact flow equation for the 
fermionic effective action IV we add an infrared cutoff 
R^ quadratic in the fermion fields 



3 an = J a p 



( F) 

Typically, the cutoff R k modifies the fermion propaga- 
tor for a momenta q 2 < k 2 such that only the fluctuations 
with q 2 > k 2 contribute effectively in the fc-dependent 
partition function. After the Legcndrc transformation 
we substract the cutoff piece in the definition of the ef- 
fective action (W = In Z) 



a/3 



(161) 



F.k 



-w + x] a i> 



k 



a/3 



(162) 



In cq. (161) j a jj are the "physical sources" and the cut- 

(F) 

off scale k is varied, with R k vanishing for k — > 0. 
The "physical effective action" is then recovered as I> = 
^F,k=o- For large k the fermionic fluctuations are sup- 
pressed and Tp,k can be evaluated perturbatively. 

An exact flow or renormalization group equation de- 
scribes the dependence of Tp,k on the scale k in terms of 
the exact fermionic propagator in presence of the cutoff 
[10] 



dkTF.kVP^] = -\dk(R k F) ) 



a 6 



r 1 (163) 

/ pa 



This can be used to extrapolate from "microphysics" at 
large k to the full (quantum-) effective action for k — > 
0. One possibility uses R^ in order to implement an 

infrared cutoff [10]. Instead of an infrared cutoff, R k F ^ 
may alternatively cut off the momenta which are close to 
the Fermi surface. For example, the temperature in the 
effective fermionic propagator may be replaced by (T™ + 
fc™) 1 /", cf. [26,13]. Approximative solutions of cq. (163) 
yield interesting results for chiral symmetry breaking in 
models with local four quark interactions [32,33,12]. Also 
the onset of instabilities in the Hubbard model has been 



discussed in this framework [22] and in a closely related 
approach [21]. 

The exact flow equation has proven to be a power- 
ful tool for renormalization group studies. One therefore 
may ask what is its correspondence in the formulation of 
the BEA. In this context we emphasize that the bosonic 
effective action is valid for arbitrary sources. It there- 
fore makes no distinction between a "physical" source j 

(F) 

and an "unphysical" cutoff term R k . In consequence, 
the addition of the cutoff term manifests itself only in the 
fc-dependence of the solution of the field equation 



QGT p -3 a p [ H k ) a[j 



(164) 



while Tb itself is independent of k. The renormalization 
group flow corresponds now to a trajectory in the space 
of solutions for fc-dependent sources! 
This may be visualized by defining 



1 



G 



a/3 



a/3 



(165) 



such that the solution Gk of the field equation corre- 
sponds to an extremum of Tj k for all fc 



or 



3,k 



dG 







(166) 



The fc-dependence of Tj k at fixed G is simply given by 
d k T lM = -\dk(R { k F) ) af} G, a (167) 

corresponding to eq. (163). This is, however, not the 
main quantity of interest. 

A central quantity is the fc-dependent location of the 
minimum of Tj k which defines the "average fermion 
propagator" Gk as a solution of eq. (166) for a fixed 
source j. The relevant exact flow equation can be ob- 
tained by taking a total fc-derivative of the identity (166), 



+ 



dG a/ 3 \q 2dG a fjdG eri 



dkG eri — 



(168) 



and reads 



dkG a , = -l(r^)~ 1 dk(R { k F) ) (169) 

A V / er/.a/J \ / afj 

The equivalence of cq.(169) with the flow of the inverse 
fermion propagator can be established by using G = 
t f } + R k F) and the identity (16), i. e. 

dkivf) =-(G- 1 d k GG- 1 ) -d k (R {F) ) 

\ J a(3 V J a/3 \ J a/3 

-l^atoF^W*, (1?0) 
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This is precisely what follows by taking two derivatives 
of eq.(163) with respect to ip, evaluated at ip = 0. 

The r.h.s. of cq. (169) involves the bosonic propa- 
gator (r^) _1 . We therefore also want to know the k- 

(2) 

dependence of T K B ' and, more generally, of the n-point 
functions corresponding to a given source j. This can be 
found by decomposing G = Gk + AG and differentiating 
r fc [AG] = Tj(G + AG) with respect to AG. As an exam- 
ple, we may investigate the /c-dependence of the inverse 
bosonic propagator 



d, 



dAG af} dAG jS J \ AG=0 

2 dGafjdG-fsdGtn \ g= q 
In a bosonic matrix notation this reads 11 



43) 



(r 



B 



(171) 



(172) 



eq. (41). More interesting is the change in the 2PI-piece 
I^ 2 ) as we will show below. 



XIII. RENORMALIZED GAP EQUATION 

In this section we derive the renormalized gap equa- 
tion (158) which has the structural form of the lowest 
order Schwinger-Dyson equation. However, using appro- 
priately renormalized sources and couplings jk and 
this equation becomes exact. 

Let us first argue that the 2PI-part (32) f( 2 ) [G] corre- 
sponds to an effective /c-dependent renormalized coupling 
Xk ■ Indeed, the expansion of T-j fc around the fc-dependent 

minimum Gk 



r 7,/t - x(rB ) [G fe ]) Q/ 3 7 iAG Q(3 AG 7 5 + 



(175) 



and similar flow equations hold for higher couplings . 
(At first sight, eq.(172) bears little resemblance with the 
fermionic flow equation. In order to establish the equiva- 
lence one needs to establish an identity similar to eq.(16) 

relating T y B ' to derivatives of Tf-) In summary, we have 
obtained a sequence of exact flow equations (167), (169), 
(172) etc. 

One could evaluate eq.(172) in the approximation (28), 
(29) where = is given by 

(T B ^)al3j5eTi = — G ae G jn G p g + G ar/ G Je G p S 

i si— I/O— l/o— 1 ^—1^—1^—1 
+ ^ae ^Sri U 0-f ~ U ari U Se U /3 7 

+ G al G /3e G 7V ~ G al G fa G je ( 1 73 ) 

and obtain 

d k (rf(G)) 

V / af: 

- G a ]d k (r[ f) ) -G P ld k (r[ f) ) 



+ G-Jdk^P) +G^dk(R k F) ) 



aS 



1^ 



+ -^G pt dt (^R k ^ G 1)a ^\ pa psG ai 

+ Xpa-a^Gpg — X pa fjjG a $ — \paasGp 1 | (174) 



{Tf[G\U l5 

(176) 

can be compared with eqs. (25), (28). This demonstrates 
that r^ 2 ) [Gk] plays the role of a running coupling and we 
define 



-G ai G [jS + G aS Gp + (f (2) [G]) Q/ 3 7< 5 



(Afe) a/ 3 7 5 - r^ 2 j 7<5 [G] 



(177) 



We emphasize that the renormalized coupling A& will typ- 
ically not describe a purely local interaction, but rather 
exhibit a nontrivial momentum dependence (form fac- 
tors). This holds even if we start from a strictly local 
microscopic interaction. 

For large values of k the renormalized coupling A& 
equals the microscopic coupling A. The lowest order 
Schwinger-Dyson equation (24) 



G af3 



-3 a (3+{Rk F) ) afj + \^ til 8G l5 (178) 



becomes always a good approximation for large enough 
k where the contribution from fluctuations is suppressed 
due to the effective cutoff R 



oo , R[ F) 



k ■ 



Indeed, for k 



oo we may write G = (—j + R 



(F) 



A)- 1 



and solve the lowest order gap equation (24) 



A 



a/3 



(-J+4" 



/ 7 5 



A = A • G (179) 



This reflects, however, only the change of the part in 
rl? which results from the term T^t = \ In det G in 



From the asymptotic behavior 12 G ~ (R^) 1 > G 1 
A • G ~ (R^)~ 2 it is apparent that the lowest order 



11 Recall that the bosonic indices correspond the fermionic 



double "indices (a/3) such that G and /,', ' ' are vectors and A ' choose R V such that a11 eigenvalues diverge for k 



-(2) 
B 

matrices. 



or are matrices, with A ■ B = \Aa_pB a p and similar for 
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Schwinger-Dyson equation becomes exact 13 in this limit. 
For large enough kg we may therefore start with the "ini- 
tial condition" Afc = A. 

We next use our definition of Afe (177) in order to write 
down an exact gap equation which is valid for arbitrary 
scales k. Indeed, defining also the "renormalized sources" 



with 



(jk)aP = 3af3 



8T 



8G. 



a/3 



[G] + -(X^afi-fsGjs (180) 



the minimum G obeys the exact renormalized Schwinger- 
Dyson equation for all k 



G a p 



-(jk)aP + {R k )a(i + ^ (Afe)Q/3 7 ^G 7( 5 



1 



(181) 



This has the generic form of the lowest order equation 
but the complexity of the problem is now encoded in the 
need to compute the renormalized source and coupling jk 
and Afc. For k —> the cutoff R y k ' vanishes. If one suc- 
ceeds to compute Afe^o and jk^o m dependence on the 
initial values A and j one ends with an exact gap equation 
relating the "physical sources" j to the full fermion prop- 
agator G. One can then use the lowest order Schwinger 
equation (21) with renormalized parameters jfc^o, Afe^ 
instead of the microscopic parameters j, A. Instead of the 
completely neglected higher order terms in a truncation 
of eq. (15) with Y a p = the approximations are now 
linked to the use of approximative values for jk and Afe . 
Even if only approximative flow equations for Afe and jk 
can be solved in practice, a large part of the contribution 
~ Y a p can be absorbed into renormalized parameters in 
this way. 

So far, the exact renormalized Schwinger-Dyson equa- 
tion (181) is only a formal organization of the field equa- 
tion (166). Its practical use depends on the ability to 
compute Afc and jk- In order to get some intuition for 
the renormalized coupling Afe we first compare it to the 
IP I- four fermion vertex Afe. Indeed, the renormalized 
four fermion coupling Afe can be related to Afe by use of 
eq. (16), i.e. 



(Afc) Q) 3 7 <5 



d A T F k 



dip a dippdip 7 dips w=0 

Fi— 1/^—1 I s^i— 1 



- r 



( r * 2) ) 



J- 1 I , ol ,a, Ga'aGppG , G 5 , 6 (182) 

(The presence of the cutoff R k does not change the 
identity (16).) In a matrix notation we write eq. (176) 
as 



vf = g + \ k = G(i + G- 1 \k) 



(183) 



1 This implies for k — > oo the lowest order expression Fj k = 
iTr In { G{R[ F) - j) } + {R[ F) - j) ■ G + i G ■ A ■ G + c(k) with 
c independent of G. 



Ga0,jS - ~ G cx^ G 0S + G al G (3j 
GafirtS ~ — G a~f G l36 + G a sGpj 

This yields the exact relation 



(184) 



(Afe) a/ 3 7 5 = Q a f3,~f& - [(1 + G X Afe) 1 Q ^a^p^f 
F^— 1 1 1 F^— 1 

= [A fc (l + S- 1 A fc )- 1 ] a/Ji74 (185) 

An expansion in Afe shows that Afe agrees with Afe only in 
lower order 

Afe = Afe — XkG~ + Afe^ _1 Afe5 _1 A/ £ — . . . (186) 

The difference between Afe and Afe is crucial in the re- 
gion where one of the composite bosons becomes mass- 

(2) 

less. Then T^> has a zero eigenvalue and we learn from 
eq. (182) that Afe diverges - as expected from the dia- 
grams with exchange of a massless boson in the partially 
bosonized language. This divergence of \k has been ob- 
served in the study of flow equations for the Hubbard 
model [20-22] or quark models models for strong interac- 
tions [32,12,33]. We emphasize that the divergence of Afe 
is inevitable if massless bosons are exchanged. However, 
the coupling Afe may remain finite even for for Afe — > oo 



Afe — (1 — XkG *) 1 Afe 



(187) 



(2) 

Indeed, a zero eigenvalue of T B suggests a finite Afe = 
I^ 2 ) in cq. (176). The important consequences for the 
running of Afe will be discussed in the next section. 

Finally, we also should acquire some intuition for the 
definition (180) of the renormalized sources jk- For all k 
we may approximate 

^7 k = n m det G + a (^k)ap-y5G a pGyS 

j ' Z a 

-\{]k) a (3G aP + l -(R k F) ) a[ G a p (188) 

Then the renormalized sources arc defined such that the 
minimum G of the approximated form (188) coincides 
with the true minimum of T-j k . The rcnormalization 
of the source term has its direct correspondence in the 
renormalization of the kinetic and mass terms in the 
fermionic language. Furthermore, using the definition 
Afe = r^ 2 ) [Gfe] the flow equation for jk is directly related 
to the one for Afe (cf. eq.(180)) 

<9fc(jfc)a/3 = -dk{^k) a pjsGjs (189) 

We will see that jk differs from j only by the effects of 
2PI-two-loop diagrams. 
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XIV. GAP IMPROVED RUNNING OF 
COUPLINGS 

Based on the results of the preceeding section we pro- 
pose in this section approximate formulae for the flow of 
Afc and j k and an approximate renormalized gap equation 
that can be used in practice. In this context we recover a 
nonpcrturbative flow equation for the renormalized four 
fermion interaction Afe that has been extensively studied 
in the past. We compare the flow equations for Afe and 
Afc and argue that the diseases of the flow of Afe are cured 
for the flow of Afe. 

Let us next first turn to the flow of the renormalized 
coupling Afe. The change of Afe reflects the higher order 
terms in T. We can write eq. (171) in the form 



With 



r^d t G = T^ t d t G + T^d t G 



(190) 



The part d k T^ t = T^ t d t G is described in eqs. (173), 
(174) and we remain with the exact flow equation <9fcAfe = 
f^dtG or 

^ = -f (») • (rf t + Afe)- 1 • d k R[ F) (191) 

Here (X k ) a/ 3^s — [G]) a p n s is interpreted as a matrix 
and 



r dct = I In dot G 



(rS 



dct) — — G n ^Gax + G n/i G 

a/3, 70 



= GaP,y8 (192) 



At this level the flow equations (189), (191) remain ex- 
act. The system of differential equations is not closed, 
however. In order to proceed we need an estimate of 
f( 3 )[G]. 

We next will derive an approximation for r*- 3 ) that 
leads to a version of the nonperturbative flow equation for 
the quartic coupling Afe which has been extensively stud- 
ied in the past [32,33,21,22]. In the same approximation 
the 2PI-character of the flow of Afe will lead to impor- 
tant improvements if the renormalized gap equation is 
exploited. As a starting point we use the definition (22) 
of Y a (i — dT/dG a p — \Xafj-ysGys and obtain the exact 
identity 



d 2 f _ , 8Y a0 

dG aP dG lS - Aaf3lS + dG lS 



= ^u(3~<& + g ^ I G a l (ry ) a p^ v 1 

= g^a/37<5 + -^pPllv^iG^G^l - GjjG^) 



(193) 



, , . , '06 ~ ^aS^ai) [ 1 B ) 
V / apa'p' OGq-' p'OG^'i/'dG^fi V / pjv'pv J 
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af3~yS 



-G al Gf35 + GasG/3^ + A Q) 3 7 5 (194) 



and eq. (173) this yields the exact expression 
8 2 f 
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+ QXpi3nv{(G a *G a g - G a lG a ^)A ap ^ u 
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d 3 T dct 
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+ {G<7<7'G p p> — G ap ' G pa ') 

d 3 r det 
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a pa' p' 



dGa> p> dGp' v < dG 1 s 
d 3 f 



dGa'p' dGp> v > dG 1 s 



^■p'v' pv 



A /r v' pv 



V / apa'p' dG a f p f dG p' i/' dG~yS ^ ' p' v' pv \ J 

(195) 

Therefore the difference between Afe and A vanishes for 
A = and f ( 3 ) = 0. 

One possible approach would be a systematic pertur- 
bative expansion for small A. Comparison with eq. (29) 
shows A ~ AG 4 as long as AG 2 remains small. In a per- 
turbative expansion the difference Afe — A is of the order 
A 2 . Indeed, the exact Schwinger-Dyson equation for T^ 3 ) 
is obtained by a further differentiation of eq. (195). We 
see that it is not closed since it involves T^ 4 \ An ap- 
proximate solution can be obtained by neglecting terms 
<~ I^ 4 ) and r^ 3 ) on the r. h. s., thus expressing I^ 3 ) in 
terms f < 2 ) , G and A. For small AG 2 one finds f < 3 ) ~ A 2 G. 

We propose here an even simpler approximation which 
is not restricted to weak coupling. For this purpose we 
take the lowest order expression (30) for Y a[} such that 



a 2 r 



dGafjdG-fS 
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~2 



[^ajpa^pSpv 



^aS pa pv)G ppG 
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Simple differentiation leads to an expression for f ( 3 ) 
8 3 f 



dGaffdGjsdGeri 
Xot5e&\f3^riv "f" X a § r ^ (7 A^j 7el /)G, 



(Aa7ecr A^3(5?^t/ A a -y^(j 

(197) 
Afe on 



After a renormalization group improvement A 
the r.h.s. this yields the flow equation for Afe 



dk(X k ) a f)-fS — -^^{X k )a~t<in{Xk) (38-qv 

-(Afe) Q ^(Afe) /37r; „}G Al ,(ri 2) ) £ - 1 crr afe (R k F) ) aT (198) 
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with (176) 

(rg } ) = -G^G" 1 +G-jG~l + (X k ) evaT (199) 

The flow equations (169) and (198) form now a closed 
system of differential equations for the running of Gk and 
Afe. The one loop character of the flow equation for Afe 
becomes obvious in the form 

dkiXkUpyS = (200) 

— ^{Xk)a-)e^{Xk) pSriv ~ (Afe) aSefj,{Xk) f} lr )v^G ^SfeG ' €r) 

Correspondingly, the flow (189) of jk is given by the "set- 
ting sun" two loop diagram. 

The exact S'D-equation (181) for the fc-dependcnt 
value of G and the exact flow equation (169) for <9feG 
are closely related - the latter follows by differentiation 
of the former using eqs. (176), (189). Instead of solv- 
ing eq. (169) one may therefore evaluate at every fc-step 
the new values of Afe and jk by the differential equations 
(200), (189). For this purpose one can insert on the r.h.s. 
of eq. (200) the lowest order 14 formula (cf. eq. (181)) 

d k G ev = -G ea d k (r[ F) ) G fin (201) 

This leads to the flow equations (159) and (160). In this 
form the close analogy to previously investigated flow 
equations for the four-fermion-coupling [20-22,32,33,38] 
is apparent. Replacing G by the "classical" average prop- 
agator (—j + yields precisely a structure similar 
to the well-known closed one loop equation for Afe. 

Actually, despite the structural similarity the one loop 
flow of Xk differs from the 2PI-flow of Xk by an additional 
term even in lowest order in Afe. Using eq. (186) we find 
in lowest order 

d k ~X k = d k X k - XkdkG^Xk (202) 
In the approximation (200) this yields 

dk(Xk)al3-,8 = —G llv dkG erj {{Xk)cx~ < en(,Xk)f}Sr 1 i' (203) 
— (Xk)a8ttj,(Xk)l3ir]u ~ (Xk )a/3ep ( Afe )jSr)v } 

and we observe that the (last) additional term as com- 
pared to dkXk (200) leads to a total antisymmetrization 
of the r.h.s. in the indices (afljS) as appropriate for a 
four-fermion vertex. Inserting on the r.h.s the lowest or- 
der relation Afe = Afe and G = (—j + R k ) reproduces 
the usual one loop flow for Afe which typically diverges 
for some finite k if the temperature is sufficiently low 15 . 



"Corrections to eq. (201) arc ~ X k G 2 d k G, cf. eq. (189) 
and could be treated iteratively. 

15 For quantitative investigations of this or similar flow equa- 
tions see [32,33] for Nambu-Jona-Lasinio type quark models 
and [20-22] for the Hubbard model. 



In the past, this divergence at finite fc has been a major 
obstacle for the exploration of the low temperature phase 
using eq. (203). 

It is possible (but not known) that the flow of Afe is 
better behaving if one replaces eq. (203) by eq. (200) 
while keeping G — (— j + Rk)^ 1 ■ Here we will advocate 
to cure the disease by a more accurate of G. Indeed, we 
want to exploit the 2PI-character of our formalism and 
propose to evaluate G by an approximation to the exact 
renormalized Schwinger-Dyson equation (181). 

In many circumstances a reasonable approximation ne- 
glects the difference between jk and j. Indeed, due to the 
2PI setting, the running of jk occurs only at two loop or- 
der as apparent from the combination of eqs. (189) and 
(159) which yields eq. (160). All one loop effects for the 
propagator-renormalization (and much more!) are con- 
tained in the "gap" ~ AG present in the implicit equation 
(181) for the approximate average propagator 

G = (-j + R k F) + X k ■ Gy 1 (204) 

Of course, the solution of the "gap equation" (204) has 
now to be evaluated for every value of k and typically 
involves additional approximations. Nevertheless, this 
procedure offers an important advantage: one can check 
at every fc-step for the possible appearance of additional 
solutions of the gap equation. In fact, in many situations 
one expects that a symmetry breaking gap A vanishes 
for large k where only a small part of the fluctuations 
is included. The flow can then produce a characteristic 
bifurcation phenomenon where below a critical scale kbf 
both solutions with A = and A / coexist. In this 
case the computation of Afe (and j k ) should follow the 
solution with lowest free energy. 

The dominant features of the term ~ A • G in eq. (204) 
can often be accounted for by a local gap as discussed in 
sect. VI, i.e. 

G = (-] + R k F) + A)" 1 (205) 

Here the local gap A a/ 3 = A a b(x)S xy involves now the 
renormalized coupling Afe . In the approximation of local 
interactions the gap equation is analogous to eq. (63) 
with A replaced by a corresponding renormalized cou- 
pling - more generally A corresponds to the local part 
of the term <~ A • G in the inverse of eq. (204), similar 
to eq. (23). Consider now the case of spontaneous sym- 
metry breaking. Once Afe has grown large enough there 
will be a solution with A^0. We define k cr such that 
for k < k cr the solution with A ^0 has the lowest free 
energy. For the running of Afe one should therefore insert 
eq. (205) into eq. (159) with A = for k > k cr and 
A^0 for k < k cr . Typically, Afe will increase towards a 
finite large value as k approaches k cr from above. There 
is no reason for a further strong increase of Afe in the 
SSB-regime for k < k cr . One therefore expects in this 
case that Afe remains finite as k — > 0, perhaps in contrast 
to an approximation which neglects A. 
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There may still remain some cases where Xk grows very 
large at some scale k s . Then, for large A^G 2 , one may 
doubt the validity of our approximations. Nevertheless, 
in many situations one may still extract useful informa- 
tion from the renormalized gap equation (181). For ex- 
ample, the value of the gap A may already settle to an 
almost fe-independent value for scales substantially above 
k s . In order to explore such a possibility one may assume 
that for k smaller than some value k the effects of the 
change of Xk and jk on the value of G are small and can 
be neglected. (This assumption is used implicitcly for all 
k for the lowest order SD-equation, whereas we use it here 
only for k < k.) We can then evaluate eq. (181) for k = 
and use the approximations j n = j k , A = X k . This 

( F} 

yields the renormalized gap equation (158) with R k 
set to zero. The combination of the gap equation with 
the renormalization flow becomes now a very powerful 
tool. A computation can proceed in two steps: First one 
solves approximately in two steps: First one solves ap- 
proximately the flow equations and computes jk , A^, . The 
second step uses the renormalized coupling and source for 
the gap equation. The validity of the approximation can 
be (partly) checked by testing if the results are indepen- 
dent of the particular choice of k. Clearly, this requires 
that has only a small influence on the value of the 
gap A. 

The renormalization group improved gap equation 
(158) should lead to a much better grasp of the physical 
situation whenever the running of couplings has an im- 
portant effect. This is certainly the case for the Hubbard 
model [20-22]. In particular, the four fermion interaction 
in the Hubbard model becomes nonlocal for k sufficiently 
below the ultraviolet cutoff (e.g. inverse lattice distance). 
As mentioned in the preceeding section, this allows for 
the solution of a nontrivial gap equation for a nonlocal su- 
perconducting order parameter. Furthermore, the error 
estimated from the dependence on k constitutes a lower 
bound for the error from the neglection of higher order 
terms ~ Y a p. More precisely, this concerns the contri- 
butions that cannot be absorbed in the renormalization 
group running. Last but not least the renormalized gap 
equation (158) offers a new perspective on the concep- 
tual use and validity of the lowest order SD-equation: 
The lowest order gap equation should be thought of us- 
ing optimally renormalized couplings and sources, such 
that residual errors are minimized. This remark is partic- 
ularly important if one is merely interested in the quali- 
tative structure of the results and less in the quantitative 
values of Afc and jk ■ The errors of neglected fluctuations 
(higher order terms) may be much smaller than naively 
estimated from the use of the unrenormalized equation! 

We close this section by a cautionary remark on the 
correct use of the renormalized gap equation. In this 
context we note that eq. (158) (or cq. (205) or similar) 
has the same ultraviolet cutoff as the leading order gap 
equation (23), i.e. the physical microscopic scale (lat- 
tice distance or similar). In particular, the scale k docs 



not act as an UV cutoff in the renormalized gap equa- 
tion, implying that fluctuations with all momenta (not 
only q 2 < k 2 ) arc effectively included. In fact, the gap 
equation often receives a dominant contribution from the 
short distance fluctuations. One therefore has to treat 
the momentum dependence of the renormalized coupling 
Xk carefully. For high "external momenta" Q the flow of 
Xk is severely suppressed if k 2 <C Q 2 . For large Q 2 one 
therefore typically finds A^ « A. On the other side, the 
running of A& may be important for Q 2 <C k 2 , modifying 
the contributions of fluctuations with small momenta in 
the gap equation. 



XV. BOSONIC RENORMALIZATION GROUP 
EQUATION 

Critical behavior is associated to an infinite correla- 
tion length for composite bosonic fields. In our formalism 
this corresponds to vanishing eigenvalues of the inverse 
bosonic propagator r^ 2) = + f ^ = T 2 dct + X k . Obvi- 
ously, the lowest order Schwinger-Dyson equation breaks 
down in the critical region since the correction 

h = (rg))" 1 - (it))" 1 
= -(rie ) t) _1 -^-(ri 2 c ) t + A fe ) _1 (206) 

grows very large. Renormalization group equations are 
the standard tool to approach the critical behavior. Even 
though formally exact the renormalized gap equation 
(181) seems not ideal for this purpose since the infrared 

(F\ 

cutoff R k acts only on the fcrmionic fluctuations. In the 
context of critical behavior it seems more appropriate to 
introduce an infrared cutoff for the composite bosons. 

Let us therefore add to the action Sj a "bosonic in- 
frared cutoff" in the form 

A k B) S=l(R ( k B) ) Mp^s (207) 
o v / a/370 

This term is quadratic in the fcrmionic bilincars in anal- 

(F) 

ogy to the piece containing the fermion cutoff R k which 
is quadratic in the fermionic fields. We will discuss 
the properties of R y k ' required for an IR cutoff below. 
Adding the term (207) amounts to a change of the ef- 
fective quartic coupling A. Defining Wk = In Zk for the 
system in presence of the cutoff and subtracting the cutoff 
again after the Legendre transform, the effective average 
action reads in the bosonic language (with j = j[G]) 

T B ,k = -W k [0,j] + \j afj G a(i - \ (r[ b) ) G a pG 7S 

(208) 

The exact flow equation for the fe-dependence of is 
now easily derived in analogy to [10] 
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dk^B.k , = -d k W k - \d k (R k B) ) G a pG jS 

| g I j ° a ^ 1& 

1 ( B) ~ ~ " " 

= \d k (R[ B) U lS {W { k 2) ) aPlS (209) 

Here we have used eq. (5). In our matrix notation we 
have = (T^ + i?^) -1 and we recover the well 

known exact flow equation for bosonic systems [10] 

d k T B , k = \tr{d k R k B) ■ (rgi + R^r 1 } (210) 

We note, however, that the trace tr acts on bosonic in- 
dices (a/3). Since the bosonic fields are bilocal it involves 
a double momentum integration. Despite the similar- 
ity in the form this constitutes an important practical 
difference as compared to the more familiar exact flow 
equation for a finite number of local bosonic fields. We 
will see below how this difficulty can be partly overcome 
for a suitable choice of R k B ^ . 

The bosonic flow equation makes no distinction be- 
tween a model with "fundamental" bosonic degrees of 
freedom and composite bosonic fields. This is a nice fea- 
ture since the conceptual difference between "fundamen- 
tal bosons" and "composite bosons" is not clear, anyhow. 
The well developed formalism of the effective average ac- 
tion and nonperturbative flow equations can therefore 
directly be used in our formulation. Furthermore, the 

(B) 

bosonic cutoff R k can be combined with a fermionic 
cutoff R k . The "initial value" of the flow for large k 

IB) 

(or k — » oo) depends on the choice of the cutoff R k 

(and possibly R y k ') and will be discussed below. 

It should be noticed that the flow actually only con- 
cerns r in eq. (41). The part Tdct remains unaffected 
since the flow equation (210) is manifestly covariant with 
respect to the transformations (34) (37). Nevertheless, 
the contribution from the fermionic fluctuation determi- 
nant is implicitely present in the flow of T 

d k f k = lt~r{d fe 4 S) (?i 2) +rW+Ri B Y 1 } (211) 
In the very simple truncation 

f k = X -G ■ A fe • G - Aj k ■ G (212) 

one has 

d k (\ k )af3~f8 = (213) 

\H d ^ \r 9 \r ^ + r det + - 
2 I K dG af idG 1 s J \g=g 

and we observe that the only contributions arise from the 
effective bosonic cubic and quartic couplings associated 
totheG-derivativesofr^ (cf. eqs. (192), (173)). There 



is, however, no need to restrict the truncation of T to the 
form (212). Then contributions to the bosonic cubic and 
quartic coupling arise also from G-derivatives of r^ 2 -* . Wc 
emphasize that the exact flow equation (210), possibly 
combined with a contribution from the ^-dependence of 
a fermionic cutoff to the running of G (cf. eq. (169)), has 
to be equivalent with the exact flow equation for partially 
bosonized versions of the fermionic model that will be 
briefly discussed in the next section. In particular, the 
partially bosonized version offers a simple starting point 
for the computation of Tb,u for large values of k. 

The simple and suggestive form of the exact flow equa- 
tion (210) can be used as a convenient starting point for 
approximations. For a practical use one has to reduce 
the complexity to a finite number of local bosonic fields 
and to truncate the corresponding bosonic effective ac- 
tion. We also have to define a suitable infrared cutoff 

( B) 

R k and to determine the initial condition for large k. 

A crucial point is the "initial condition" for the flow for 
k = A. (Instead of k = A one may also use k — > oo.) One 

( B) 

would like to choose R A such that Tb,a can be com- 
puted reliably and takes a simple form. A particularly 
interesting choice is (in a matrix notation) 

R (B) = -A (214) 

In this case the bosonic cutoff precisely cancels the inter- 
action term such that Wa is given by a free theory. Since 
the cutoff is substracted again in the definition (208) of 
T^A one has the exact solution 

r B ,A - Ti (215) 

At the scale A the lowest order SD-equation becomes 
therefore exact! 

The missing fluctuations are next included by "switch- 
ing on" the coupling in the higher order fluctuation effects 
as k is lowered from A to zero. For RIZ'q the full BEA 
- including all fluctuations and all interaction effects - is 
recovered from the solution T B ^ . 

( B) 

We next have to define a suitable cutoff R k for in- 
termediate k. For this purpose it is useful to interprete 
G a fi as an infinite number of bosonic fields labeled by 
the relative coordinate z = y — x and to perform a 
Fourier transform with respect to the "center coordinate" 
s = (y + x)/2. In this basis one has 

G a[j - G ab (Q, z)= f e-^ s G ab (s, z) (216) 

and the local gap discussed in sect. VI corresponds to 
the field labeled by z — 0,G a b(s,z) — g a b(s)S(z). For 
local interactions it is sufficient to have the cutoff acting 
only on g a b(s) and we preserve translation invariance by 

(R k B) )ab,cd{Q,Z,Q',z') = 

j j e iQz e -iQ'z' (jR ( s )) a6cd ( S) ZjS >, z >) 

J Z J z' 

= (R (B) )a bcd (Q)S(Q - Q')5(z)6(z') (217) 
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with (R A )abcd{Q) = —^abcd- We note that for this cut- 
off the trace in eq. (210) reduces to a single momentum 
integral (plus summation over internal indices). 

For the minimal set of spinors, X abc d = — Ae a bcd) the 
initial effective action T\ depends 16 only on six complex 
bosonic fields g a b(s). By virtue of the discrete symmetry 
$ (sect. V) this can be reduced to four real and one 
complex fields. With respect to the SO(4)=SU(2) s x 
SU(2) C symmetry transformations these fields transform 
as (3, 1) + (1, 3). The triplet m with respect to the spin 
transformations 17 

m( s ) = (0( s )t?/>(s)) , 
mi = #32 + 541 , m 2 = -i(032 - 54i) , 
m 3 = #31 - 542 (218) 

and the triplet c with respect to the "charge" transfor- 
mations SU(2) C 

ci = -(512 - 534) , c 2 = -2(512 + 534) , 

c 3 = @i/>) = .gal + 542 (219) 

span the irreducible representations. The third gen- 
erator of SU(2) C corresponds to electric charge. The 
charged boson corresponds to the electron pair 512 = 
(ipiT2ip) = — {c\ — ic2)j2 with an associated hole pair 
534 = (i>iT2ip) = (ci + ic 2 )/2. Labeling by s,t the irre- 
ducible bosonic representations we write 

R^\Q 2 )=Kr s (^^)s st (220) 

where r s (k = A) = 1 and A s are determined such that 
— — A. We also require that r s vanishes fast for 
k 2 <C Q 2 such that the fluctuations with high Q 2 are 
already effectively included in Tg^. Here Q s is the de- 
viation of the momentum Q from the momentum Q s for 
which the free energy of the boson s is minimal 18 . 

As a possible truncation for practical computations we 
propose to keep two sets of local bosonic fields. One set 
is selected from the g a b{s) whereas the other set corre- 
sponds to local gaps A ab (s) according to the ansatz 

G- b \s,z) = -Uk)ab(z)S(s) + A ab (s)d(z) (221) 

This allows us to express T^et as a functional of A, and 
similar for in the flow equation (211) for IV For Tk 
we make the ansatz 

f fe = T' k [g, A] - \[{j k ) a p - j a p]G aP (222) 



and derive the flow equation for j k from the first G- 
derivative of eq.(211), evaluated at s = , z ^ 0. The 
remaining flow equation 19 for T' k 

dk T' k = ltr{d k R k B) (T' k (2) + lf e >[A] + R^)- 1 } (223) 

obtains by restricting eq. (211) to fields depending only 
on s. Corresponding to the choice of cutoff (217) the 
trace tr involves now only one momentum integration 
and summation over the internal bosonic indices (ab) . At 
this stage, the inversion of the matrix of second functional 
derivatives has still to be done in a larger space, however. 

We have now arrived at a flow equation for r' fc which is 
already rather similar to the well studied flow equation 
for a finite number of bosonic fields. The truncations 
developed in this context, e.g. a derivative expansion 
[34] can be applied. The only unusual contribution is 

the part ~ T^ t [A] which accounts for the fermion loops 
in presence of the gap. 20 In order to relate 5 and A 
one minimizes the free energy Tj = F/T, noting that the 
relation 

(jk)ai3G aP = / A ab (s)g ab (s) + const (224) 

J s 

is independent of k and therefore 

^±=T' k [g,A]- 1 - J A ab (s)g ab (s)-^Trln(-j k + A) + c 

(225) 

In particular, in an approximation where V is indepen- 
dent of A the variation eq. (225) with respect to 5 yields 
a simple relation analogous to eq. (76) 

dr' 

A ab (s) = j—f- (226) 
dgab(s) 

This expresses A as a functional of 5 and allows us to 
study the flow equation (223) for V being a functional 

only of the local bosonic fields g a b(s)- Nevertheless, the 

(2) 

term r^ et still involves the source jk according to eqs. 
(192), (221). _ 

The situation becomes particularly simple for the 
study of critical phenomena in a context where all 
fcrmions acquire a nonzero mass (gap). Then for low 

enough k we may neglect the fermion fluctuations alto- 

( 2) 

gether and omit the term in the flow equation (223). 



16 This property will not be preserved for F fe , k < A, due to 
the contribution from in the flow equation (211). 

17 We omit the tilde in the notations for the spinors here. 

18 If the free energy at a given scale k is minimal for a ho- 
mogenous field g s one has Q s = Q. For a preferred antiferro- 
magnet one would rather have Q a = Q — Q a ■ 



19 Note that the second term in eq. (222) does not contribute 
toff. 

20 Note that in absence of a fermionic cutoff R k F ^ the 
quadratic part of the fermion fluctuations is already included 
in r^A (215). On the other hand, for -R A F ' — * 00 the mini- 
mum G is in a region where fermion fluctuations are effectively 
absent for k = A. 
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What remains is the well studied flow equation for a 
bosonic system that has already given a very satisfactory 
description of critical phenomena [10,34]. This demon- 
strates how critical phenomena can be understood in the 
context of the BEA. If massless fermions are present at 
the phase transition (as, e.g., the antiferromagnetic tran- 
sition in the Hubbard model [13]), the universal critical 
exponents will be affected by the fermionic fluctuations 

(2) 

and Tfe t cannot be neglected even near the phase tran- 
sition. 

Let us summarize sects. XII- XV from the viewpoint of 
practical applications. Employing suitable fermionic and 
bosonic cutoffs and approximate truncations we have re- 
covered the "standard" well investigated approximations 
for non-perturbative flow equations. This holds both for 
the flow of fermionic and bosonic systems - both are de- 
rived from a common exact renormalization group equa- 
tion for the BEA. Furthermore, our 2PI setting points to 
extensions of the simplest versions of these "standard" 
flow equations. For the fermionic flow equation the lead- 
ing effect of the effective bosons can be incorporated via 
the inclusion of the gap in cq. (200), using at every k 
an approximative solution of eq. (202). For the bosonic 
flow equation the effect of the fermionic fluctuations is 
incorporated via the inclusion of T^JA] in cq. (223), 
with A determined by eq. (226). However, even for 
jk = 3 an d using eqs. (192) (205) for T y d J t , the fermionic 
contributions still appear in a rather complicated form. 
At the present stage the bosonic renormalization group 
equation may therefore be mainly considered as demon- 
stration how the well known bosonic critical behaviour 
arises in a fermionic model. For practical investigations 
of the phase transition the functional renormalization in 
the partially bosonized formulation (next section) may 
be easier to handle. 



XVI. FUNCTIONAL RENORMALIZATION FOR 
PARTIAL BOSONIZATION 

In the formulation of partial bosonization (PB) the si- 
multaneous implementation of cutoffs for fermions and 
bosons is particularly simple. One adds to the action 
in eq. (145) a piece ~ -ft! quadratic in the fermions 

ip and another one ~ W k ' quadratic in the bosons g. 
Within PB the derivation of the exact renormalization 
group equation is straightforward [10] 

d k T k = ^Str ^d k R k (T^ + J R fe )" 1 | (227) 

Here stands for TpB,k as defined by the Legendrc 
transform in presence of the cutoff. Again, the cutoff 
term is subtracted in generalization of cq. (162) (or eq. 
(208)). The supertrace Str contains a minus sign for the 
fermions (cf. eq. (163)) and involves now only one mo- 
mentum integration according to a one loop expression 



for local fields. The matrix acts in the combined 
space for fermions and bosons, just as R k . Truncations to 
cq. (227) have been successfully investigated for Nambu- 
Jona-Lasino type models [12], the Gross-Neveu model 
[38] and the Hubbard model [13]. 

We may compare eq. (227) with the fermionic flow 
equation (sects. XII- XIV) and the bosonic flow equation 
(sect. XV). For the fermionic flow equation we have ar- 
gued that a major improvement for the treatment of the 
low temperature phase may arise from the inclusion of the 
gap (cf. eqs. (200) (204)) and the concentration on the 
coupling Afc instead of the four fermion vertex X k - Both 
features are naturally implemented by eq. (227). The 
gap is present by the Yukawa coupling ~ ipipg. Typically, 
the flow follows the effective average potential U k {g)- In 
case of spontaneous symmetry breaking a minimum of U k 
for g ^ appears at some value of k. Just as discussed 
for the gap improved renormalization group equation in 
sect. XIV we can easily follow the minimum of Tj^fb 
in the regime with spontaneous symmetry breaking and 
thus account for a nonvanishing gap. 

The flow in the bosonic sector typically monitors the 
inverse bosonic propagator P for a few local fields gi=g ab . 
Let us assume a basis where P is diagonal and denote the 
eigenvalues by Pi (q) (in momentum space) . Furthermore, 
one may compute the flow of the corresponding Yukawa 
couplings 

A/ = - / H iab{q, q')Ml)Ml')9i{-<l - <?') (228) 
Jq-q' 

In the simplest approximation one may keep a fixed mo- 
mentum dependence such that only the constants h iab 
depend on k 

H tab (q,q') = h lab F tab (q,q r ) (229) 

whereas the momentum dependence in Fi ab (q,q / ) de- 
scribes a given channel, like antiferromagnetism in sect. 
VIII. We can now reconstruct the part of the effective 
four fermion coupling that is induced by the exchange of 
the bosons. For this purpose one solves the field equation 
for gi as a functional of tp 

+Pr\q) f H iab (q-q',q')ip a (q-q')ip b (q') 

Jq> 

+ ■■■ (230) 

Inserting this solution into T^fb and extracting the term 
involving four fermion fields yields the contribution to the 
four fermion vertex. For vanishing (g^ = 0, Ji = this 
reads 

z i Jq,q',q" 
H iab (q-q',q')H icd (-q-q",q") 

■Mq - q')Mq')M-q - q")Mq") (23i) 
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We recognize the close analogy to the "boson exchange 
contribution" ~ (T^) -1 in eq. (182) and recall that A^ 

is related to by eq. (183). The information contained 
in Pi(q) and H ia b(q,q') contains an essential part of the 
information related to A&! 

An important difference of principle, remains however. 
This arises from the restriction to a finite number of local 
boson fields (x) instead of the bilocal bosonic quanti- 
ties G a b(x,y). In consequence, the most general form of 
the four fermion vertex cannot be obtained by eq. (231). 
The contribution (231) is a sum of terms that factorize 
in the form f(q 1 , <7 2 )/(<73, ~9i ~92 ~93)M9i +92) whereas 
the four fermion vertex Xk(qi, 92, 93) can be an arbitrary 
function of three independent momenta. Of course, the 
difference can be accounted for by the remaining piece of 
the four fermion interaction which is 1PI in the fermion- 
boson-model. Despite this limitation, a truncation of the 
effective four fermion vertex to the sum of factorizable 
terms (231) is sufficient for many purposes. Simple ap- 
proximations can still describe a rich momentum struc- 
ture of Afc by a finite number of parameters consisting of 
h ia b as well as masses and wave function renormalizations 
for the momentum dependence in Pi(q). Obviously, in- 
creasing the number of boson fields gt(x) retained in PB 
improves the resolution of the momentum dependence of 
Afe. 

At this point we may recall the redundancy within the 
fcrmion-boson model: a given piece in A^ with the struc- 
ture (231) can both be obtained from the proper vertex 
(1PI in the fermion-boson-model) and from the exchange 
of the bosons Redundancies of this type should be 
avoided in order to keep the number of running couplings 
manageable. After all, the split into different pieces con- 
tains no relevant information. It is indeed possible to 
shift all contributions of the type (231) uniquely to the 
exchange of bosons at every scale k. For this purpose the 
corresponding contribution from the proper vertex which 
is generated by the flow can be transferred to the bosonic 
sector by the method of "rebosonization" [11]. This uses 
a fc-dependent definition of the bosonic variables. Re- 
bosonization is a crucial ingredient for truncations that 
omit the proper four fermion vertex! 

Even for simple truncations the functional renormal- 
ization in PB is a powerful method to include the bosonic 
fluctuations beyond mean field theory. This has been 
tested by studies of the mean field ambiguity [8,13]: For 
phase transitions the dependence of the critical tempera- 
ture on the choice of partial bosonization (S a bcd) is often 
very substantial in mean field theory. It turns out to 
be greatly reduced once the effect of the bosonic fluctua- 
tions is included via a solution of the flow equations. For 
the two dimensional Hubbard model with next neighbour 
interactions the partially bosonized flow equations yield 
already a rather satisfactory picture at small doping [13]. 
One finds a pseudo-critical temperature T pc which is con- 
siderably below the value of fig. 1 due to the effect of the 
bosonic fluctuations. Furthermore, T pc is not the critical 



temperature for the transition to antiferromagnetic or- 
der. For a temperature range T c < T < T pc the system 
is in the symmetric phase (vanishing antifcrromagetic 
order parameter) and the correlation length increases 
- ^ 1 exp(-6T c /T),6 w 20.7. This describes the uni- 
versal critical behaviour of the classical (3(3)-Heisenberg 
model [35] and reflects that the evolution for small k is 
dominated by the Goldstone bosons. For k <C T < T pc 
the flow in the fermion-boson-model reduces to a simple 
flow of a bosonic linear cr-model [10], [31]. Finally, for 
T ~ T c the correlation length reaches the macroscopic 
sample size L w 1cm. For T < T c long range antiferro- 
magnetic order occurs and the temperature dependence 
of the order parameter has been computed quantitatively. 
The order becomes possible due to the absence of fluctu- 
ations with wave length larger than L in a finite sample. 
The critical temperature shows a weak logarithmic de- 
pendence on the sample size and vanishes for L — > 00, 
in accordance with the Mermin- Wagner theorem [36]. 
All this demonstrates the practical viability of functional 
rcnormalization in a partially bosonized setting. Similar 
achievements for the flow equations discussed in sects. 
XII- XV have still to be demonstrated. 

In the present paper we have developed the formal im- 
plementation of the functional renormalization group for 
the BEA. Concrete examples have to be investigated in 
order to assert the practical use of the proposed approx- 
imate flow equations. Nevertheless, we have already es- 
tablished a rather close contact with the functional renor- 
malization group in an approach with partial bosoniza- 
tion. The success of the latter is a strong argument in 
favor of the conclusion that the main features of renor- 
malization can indeed be obtained from rather simple 
approximations to the exact flow equations. 

XVII. CONCLUSIONS 

We have discussed a bosonic effective action (BEA) for 
interacting fermion systems. While the general form de- 
pends on bilocal fields or, equivalcntly, infinitely many 
local fields, the structure becomes much simpler for 
suitable approximations. The lowest order of a sys- 
tematic loop expansion is equivalent to the lowest or- 
der Schwinger-Dyson (SD) equation. As an important 
achievement, the BEA offers a simple and direct way for 
the computation of the free energy corresponding to dif- 
ferent solutions of the SD-cquation. This is crucial in 
order to establish the ground state in such a situation. 
Furthermore, for local four fermion interactions it is suf- 
ficient to consider local gaps in the lowest order of the 
loop expansion. Therefore the BEA only needs a finite 
number of local fields in this approximation. Already at 
this level useful computations can be performed as we 
have demonstrated for the phase diagram for the Hub- 
bard model with next-neighbour interactions in two di- 
mensions. We emphasize that the one loop BEA goes far 
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beyond the perturbative loop result! In particular, it can 
describe spontaneous symmetry breaking. 

Going beyond the leading order gets more involved. 
For the BEA it is a priori not obvious how a restriction 
to a finite number of bosonic fields can be implemented 
in a systematic way. We propose here to use a truncation 
of an exact renormalization group equation. For this pur- 
pose we generalize the exact renormalization group equa- 
tion for the (1PI) effective average action [10] to a 2PI 
-formulation. It describes the dependence of the BEA on 
a characteristic length scale The price to pay for two 
particle irrcducibility is the appearance of bilocal fields 
Gap. As long as only a fermionic cutoff is introduced for 
the description of the renormalization group running and 
the truncation is kept to low orders of G one can cope 
with the problem of bilocality. We propose to use a sim- 
ple renormalized gap equation (158). The renormalized 
couplings appearing in this equation obey rather simple 
approximative flow equations (159) which can avoid the 
previously observed divergence of the couplings. 

A reliable description of the flow for critical phenom- 
ena or other situations with a large bosonic correlation 
length calls for the additional introduction of a bosonic 
infrared cutoff. By a suitable choice of the cutoff Rk we 
can reduce the right hand side of the bosonic flow equa- 
tion to a sum of single one-loop momentum integrals over 
the propagators of a finite number of bosonic fields. The 
infinite number of bosonic fields contained in G a p ap- 
pears now only in the relation between the propagators 
and the functional derivatives of the BEA. We propose a 
truncation that reduces this problem to a finite number 
of bosonic degrees of freedom as well. In summary, the 
flow equation interpolates from the solution of the lowest 
order gap equation at short distances (k = A) to the full 
2PI quantum effective action as the cutoff is removed for 
k — > 0. The practical use of this particular approach has, 
however, to be demonstrated in the future. 

The discussion of the functional renormalization in the 
2PI-contcxt of the BEA is also useful for a direct com- 
parison with formulations based on partial bosonization. 
We find that already simple truncations for the partially 
bosonizcd functional renormalization flow reproduce the 
main improvements of the 2PI-flow in a natural way. 
The simple formulation of the partially bosonized flow 
in terms of a finite number of local field appears to be an 
important advantage for practical computations. 

At this point we may compare the merits and disad- 
vantages of the BEA to the other competing method to 
go beyond mean field theory, namely partial bosonization 
(PB). Let us first look at the criteria formulated in the 
introduction, (a) Both BEA and PB permit a reason- 
able simple computation of the free energy (b) Whereas 
BEA is one loop exact in lowest order by construction, 
this property is realized for PB only after inclusion of 
the bosonic fluctuations, as demonstrated in a RG ap- 
proach [8]. (c) Both BEA and PB can include the ef- 
fective bosonic fluctuations in a more or less systematic 
way by the use of renormalization group equations. The 



BEA offers the additional advantage of a simple estimate 
of the importance of the bosonic fluctuations from a com- 
putation of (parts of) the three loop effective action (see 
sect. VII). On the other hand, the flow equations for PB 
are both conceptually and practically much simpler. 

The BEA has the important merit of avoiding redun- 
dancy in the description from the outset. For PB this 
redundancy leads in lowest order to the problem of the 
Fierz ambiguity. It can, however, be (partially) absorbed 
by rebosonization [11], [8] in the flow equations. The 
Fierz ambiguity can even be turned into an advantage 
by providing an error estimate for approximations. An- 
other asset of the BEA is the high degree of symmetry 
which allows one to connect different channels. Indeed, 
the BEA is a very flexible approach since the various 
channels where condensates or other interesting physics 
may occur are all described at once. This feature is very 
helpful for the exact mapping between the repulsive and 
attractive Hubbard model presented in this paper. In 
contrast, PB needs a selection of the interesting channels 
from the beginning of its formulation, thereby introduc- 
ing a certain rigidity and perhaps sometimes a certain 
bias. On the other hand, this restriction of PB to a few 
degrees of freedom is an important advantage for practi- 
cal computations. Comparing this last aspect with BEA 
we find no practical problem in lowest order. In higher 
orders, however, an effective reduction to a finite num- 
ber of bosons requires some thought as discussed in sects. 
XII-XV. 

Finally, PB has the advantage that both bosons and 
fermions appear explicitely. This permits a simple de- 
scription of universality classes for (near) critical be- 
havior for which only the fermionic and bosonic low- 
mass fluctuations are relevant. The use of nonpertur- 
bative flow equations for fermionic systems with com- 
posite bosons has been well tested in this context and 
has reached a highly sophisticated level. For BEA the 
fermions appear not explicitely anymore. Nevertheless, 
the fermion fluctuations are included and the running of 
the "fermionic part of the effective action" is reflected by 
the running of jk (sects. XII, XIII). It may be advanta- 
geous to cast the fermionic fluctuation effects into a more 
explicit form. This, as well as a practical demonstration 
of the use of the flow equation remains to be done in the 
setting of BEA. 

Before a concrete computation of the flow of the BEA 
for a simple fermionic system has been performed it is 
perhaps premature to judge the relative merits of the two 
approaches for systems with a long correlation length. It 
seems well conceivable to us that both the bosonic ef- 
fective action and partial bosonization should be used si- 
multaneously in order to assert how reliable are proposed 
solutions for strongly correlated fermionic systems. Such 
an approach can use the presented detailed relations be- 
tween the quantities computed by both methods. 
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APPENDIX A: GROSS-NEVEU-MODEL 

The action for the Gross-Neveu model (GN-modcl) [23] 
can be written in the form 



S = Jdrj d D x[^d^ + ^(?V) 2 } (A.l) 



For D < 3 and (3 — > oo this action is invariant under 
the rotation ("Lorentz") group SO(D + 1) and there- 
fore describes relativistic fermions. This holds provided 
the matrices 7 M are appropriately chosen, where /i = 0, 
1 . . . D runs over d = D + 1 values. For D = 2 we take 
7° = T3, 7 1 = n, 7 2 = T2- An appropriate choice for 
D = 3 is 7 = 1, 7 1 = in, 7 2 = «T2, 7 3 = 1x3 and we 
note that -0 accounts only for one chirality, i. e. ip=ipL, 
ip=ip L . The 50 (4)-invariance of this somewhat unusual 
"Chiral Gross-Neveu model" can be seen by noting the 
"Fierz identity" for two component spinors 

«¥) 2 = -\{^){^) = ^v^xv^vo (a.2) 

Finally, for D = 1 we may take 7 = r 2 , 7 1 = n- 

The continuous GN- model (A.l) needs some type of 
an ultraviolet regularization. A cutoff preserving the 
SO(d) symmetry can easily be implemented in momen- 
tum space by restricting q^q^ < A 2 . Another approach 
provides for a suitable smooth momentum dependent 
cutoff by modifying the fermion kinetic term in the ac- 
tion (A.3). As an alternative, we may discretize the 
model on a cubic lattice with lattice distance a such 
that ip(x,r) = a~ D ^ 2 tp(mi, 7712, . . . mo, t), with integers 
(mi,m2,...mn) denoting the lattice sites. Depending 
on the precise implementation of the discretized deriva- 
tive the continuum limit a — > of the lattice GN-modcl 
may correspond to a continuous GN-model with several 
species of fermions ( "fermion doubler problem" ) . 

In the following we concentrate on D = 2. On a cubic 
lattice with sites (m, n) the action of the lattice-Gross- 
Neveu model (LGN-model) 

S = E {^n + \u{^) 2 ) , 

(m,n) 

Uin = i?7 M ^V (A.3) 
is formulated in terms of the lattice derivatives 

diip(m, n) = (^ip(m + 1, n) — tp(m — 1, n)j 

d2ip(m, n) = — (i>(m, n + 1) — ip(m, n — l)j (A. 4) 

with da = d T and 7 = T3, 7 1 ' 2 = n,2- The coupling 
strength is related to the one of the continuum formula- 
tion by U = G/a 2 . 



Without changing the structure of the interaction term 
~ U we may bring the discretized action (A.3) of the 
GN model closer to the one for the Hubbard model (55) 
by appropriate redefinitions of the spinor fields. This 
transformation is an example for the extended symmetry 
transformations which leave Tb invariant, as discussed in 
sect. IV. With 

urn 

tp(m,n,T) = -ir 3 cxp(-— n) 

iirm , 
exp(-^-r 2 )V' (m,n,r) 

V>(to, n, t) = ip '(m, n, t) 

iirm inn 

cx p( — y T2) exp (~ Tl ) ( A - 5 ) 

the kinetic term equals the one for the Hubbard model 
up to an important factor (—1)™, i. e. 

4m = ip '{m, n, T)d T ip'{m, n, t) - tip '(m, n, r) 

{ (-1)™ (i>'{m + 1, n, t) + y/(m - 1, n, t)) 

+V>'(m,n+l,r) + V'(m,n-l,r)| (A.6) 

We identify t = — l/(2a) and note that a change of sign 
of t can be achieved by further multiplying both ip'(m, n) 
and ip '(m, n) by a factor (-l) m+ ™. 

The alternating factor (—1)™ multiplying the next 
neighbour interaction in the 1-direction cannot be ab- 
sorbed by further local transformations of the fields. (It 
may be shifted to a factor (— \) rn in the 2-dircction.) It 
expresses the difference in the structure of the Fermi sur- 
face between the LGN- and the Hubbard model. We 
may call the model with the kinetic term (A.6) the "lay- 
ered Hubbard model" (LH-model) since the next neigh- 
bour interaction in the 1-direction switches sign between 
different layers in n. The hopping of the electrons in 
the 1-direction differs between the layers with n even or 
odd. In physical terms, such a situation may be realized 
by an appropriate atomic lattice structure. Results on 
phase transitions in the Gross-Neveu model [37] , [38] can 
be carried over to the LH- model. 21 Similarly, by use of 
the inverse map (A. 5) the Hubbard model (55) can be 
mapped onto a "layered Gross-Neveu" model where the 
factor (—1)™ now multiplies the derivative term ~ T\d\. 

The exact map between the LGN and the LH mod- 
els opens new perspectives for the mutual understanding 
of both models. For example, the (pseudo-)scalar order 
parameter in the LGN-model appears as the third com- 
ponent of an antiferromagnetic order parameter in the 
LH model 



21 We expect the LGN-model to belong to the universality 
class with four independent two component spinors in a con- 
tinuum formulation (N = 4). This is connected to "fermion 
doubling" on the lattice. 
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iip(m,nT)ip(m,n,T) = 
{-l) m+n ^P '(m, n, T) T3 iP'(m, n, r) = a 3 (A.7) 

The LH model is invariant under independent global spin 
rotations ip' — > exp(iQ?r)'0', ip ' ^> ij) 'exp(— iaf). By 
the mapping (A. 5) these symmetries appear as (to, n)- 
dependent transformtations in the LGN-modcl. In par- 
ticular, the other two components of the antiferromag- 
netic spin vector in the LH model appear in the LGN 
model as 

ai = {-i) m+n ^ '(to, n, T)nip'(m, n, r) 

= (—l) n tp(m,n,T)T2'ip(m,n,T) 
a 2 = {-l) m+n ^ '(to, n, t)t 2 /(to, n, r) 

= (-l) m+1; 0(TO, n,T)n-0(TO, n,r) 

(A.8) 

From the symmetry of the LH model we infer that the 
free energy in the lattice Gross-Neveu model is degener- 
ate in the direction of the vector (di, 0,2,0,3). This would 
probably not have been suspected from a direct inspec- 
tion of the action (A. 3). In particular, we predict for the 
ordered low temperature state of the lattice GN-modcl 
the existence of two massless Goldstonc bosons. 
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Note added 

The formalism presented here has been generalized to 
higher multi-fcrmion interactions in [18]. In this paper 
chiral symmetry breaking in strong interactions (QCD) 
has been studied for a six fermion interaction based on 
instantons. The functional renormalization group in the 
partially bosonized version has been applied to an inves- 
tigation of the antiferromagnetic order parameter in the 
Hubbard model in [13]. The material presented in sects. 
IX- XI and XVI was not included in a first version of this 
article. 
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